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Abstract 

We define a geometric flow that is designed to change surfaces of cylindrical 
type spanning two disjoint boundary curves into solutions of the Douglas-Plateau 
problem of finding minimal surfaces with given boundary curves. We prove that 
also in this new setting and for arbitrary initial data, solutions of the Teichmuller 
harmonic map flow exist for all times. Furthermore, for solutions for which a three- 
point-condition does not degenerate as t oo, we show convergence along a se¬ 
quence ti —>■ oo to a critical point of the area given either by a minimal cylinder or 
by two minimal discs spanning the given boundary curves. 


1 Introduction 


Tcichmiiller harmonic map flow, introduced in the joint work [12] with Peter Topping for 
closed surfaces, is a geometric flow that is designed to change parametrised surfaces into 
critical points of the area. Indeed, for closed surfaces in a non-positively curved target 
manifold, the flow always succeeds in changing, or more generally decomposing, the 
initial surface into (a union of) branched minimal immersions through globally defined 
smooth solutions [15]. 

Here we take a first step to generalize this approach to the problem of flowing surfaces 
with boundaries to a solution of the Douglas-Plateau problem of finding a minimal surface 
spanning given boundary curves. Namely, given two disjoint, closed C 3 Jordan curves 
r± in Euclidean space we investigate how to flow a surface of cylindrical type in order 
to find a minimal surface spanning the two boundary curves. 

As in [12] this flow will be constructed as a gradient flow of the Dirichlet energy 



considered as a function of two variables: a map u : Cq —> K” parametrising the evolving 
surface over a fixed domain, here the cylinder Cq = [—1,1] x S 1 , and a Riemannian 
metric g on the domain. 

We remark that if a pair ( u , g) is a critical point of E then it is also a critical point 
of the area, to be more precise either a constant map or a (possibly branched) minimal 
immersion [7]. A key idea of Teichmuller harmonic map flow is to consider E on the set 
of equivalence classes of maps and metrics modulo the symmetries of E , compare [12]. 
Thus we identify 
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• ( u , g) ~ (u, A • g) for all functions A : Co —> R + due to the conformal invariance of 
E, and 

• (u, g) ~ (u o /, f*g) for diffeomorphisms / : Co —> Co that are homotopic to the 
identity. 


As in [12] we shall then define Teichmiiller harmonic map flow from cylinders as an L 2 
gradient flow on the resulting set of equivalence classes 

A := {(it, < 7 ) : u : Co —> K" so that it|{±i}xSi parametrises r±, g a metric on Co}/ ~ . 

One important point to be understood in order to truly define such a flow is how to define 
an L 2 -metric on A and the closely related question of how to best represent a curve in 
A through pairs of maps and metrics. 

In the definition of Teichmiiller harmonic map flow on closed surfaces in [12] we chose a 
canonical representative by asking that g has constant curvature K g = 1,0, —1 (depend¬ 
ing on the genus) and that || 9 t p||i 2 is minimal. For the resulting L 2 -gradient flow this 
means that the symmetries are used to maximally simplify the evolution equation for the 
metric component in the sense that g only moves by the part of the gradient of E that 
is orthogonal to the action of the symmetries. At the same time, the map component 
evolves with the full gradient, i.e. the tension. As a result, cf. [13], for closed surfaces 
the evolution of the metric turns out to be well controlled as long as inj (M,g(t)) 0, 

while the map component shows a similar behaviour as the solutions of the corresponding 
flow for fixed metrics, i.e. the harmonic map heat flow of Eells-Sampson, which is well 
understood for closed domain surfaces and maps into general target manifolds. 

In the present setting of flowing surfaces with boundary the situation is somewhat differ¬ 
ent, mainly because our boundary condition is not of Dirichlet-type, but only of Plateau- 
type, i.e. prescribing the boundary values only up to reparametrisation. As such, even 
for a fixed metric go, one cannot expect strong regularity results for the gradient flow of 
u 1 ^ E(go , it) unless one imposes a three-point-condition for u\gc 0 - 

For maps u parametrised over the disc (Di(0), g euc i) such a gradient flow of maps was 
introduced and studied by Chang and Liu in [1, 2, 3] who considered both maps into 
Euclidean space and into Riemannian manifolds. The more general case of flowing to 
discs of prescribed mean curvature (and prescribed Plateau-boundary condition) has been 
considered more recently by Duzaar and Scheven [ 6 ]. They show that an isoperimetric 
condition on the prescribed mean curvature ensures the existence of global weak solutions 
and that these solutions subconverge to a disc with the prescribed mean curvature. 
In both cases, the flows are given as equations for only a map component u. This is 
consistent with our approach as the special structure of the disc makes it unnecessary to 
also evolve a metric on the domain; namely, the moduli space of the disc consists of only 
one point and one can furthermore pull-back any map by a suitable Mobius transform to 
obtain a map that obeys a three-point-condition. Since Mobiustransforms do not change 
the conformal structure this means that one can replace ( u,g euc i ) by a representative 
of the same point of A whose map component satisfies a three-point-condition without 
having to adapt the metric component at all. 

These special features of maps and metrics on the disc are not present for any other 
surface with boundary, though in case of the cylinder the moduli space has a very simple 
structure as it is one dimensional. But even in this case, the group of conformal diffeo¬ 
morphisms from (Co, g) to itself is not sufficiently large to impose a three-point-condition 
for the map component without having to adjust the metric suitably. As some kind of 
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restriction on how u\qc 0 parametrises the boundary curves r± is needed to obtain a flow 
that admits global solutions, we shall thus use the symmetries in a slighly different way 
than in the case of closed surfaces. Namely, we use only most, but not all, symmetries 
to ensure that the evolution of the metric is regular, while also setting aside a number of 
degrees of freedom (3 per boundary curve) to prevent a formation of singularities of the 
map at the boundary by imposing a three-point-condition. 

We remark that the Douglas-Plateau problem has been considered by many authors and 
we refer to the books [5], [4], [9], [17] and the references therein for an overview of existing 
results. What we would like to point out is the well known fact that while one can in 
general not prescribe the topological type of a minimal surface, one always obtains a 
minimal surface that is parametrised either over the original domain, for us the cylinder, 
or over a surface of a simpler topological type, in the present situation two discs. 

The paper is organised as follows. To begin with, we discuss how to best represent 
curves in the set of equivalence classes A and consequently give the precise definition of 
the flow. We then state our main results which guarantee the existence of global weak 
solutions for arbitrary initial data, see Theorem 2.6, as well subconvergence to either a 
minimal cylinder or to two minimal discs spanning the given boundary curves, at least 
for solutions for which the three-point-condition does not degenerate, see Theorem 2.7. 
The rest of the paper is then dedicated to the proof of these results. In section 3 we prove 
short-time existence based on a time-discretisation scheme and derive a priori estimates 
on the map and metric component which are crucial for both the proof of existence and 
of asymptotic convergence. This asymptotic analysis is carried out in section 5 but before 
that, in section 4, we establish that solutions exist for all times. 


2 Definition of the flow 

2.1 Representing a curve in A: Admissible variations 

As preparation for the definition of Teichmiiller harmonic map flow on cylinders we 
discuss ways of representing curves in the set of equivalence classes A through suitably 
chosen pairs of maps and metrics. We do not claim that our choice is canonical but 
rather that it is designed for the purpose of obtaining a gradient flow of energy that 
admits global regular solution. 

To begin with, we need to identify a suitable representative of a conformal class c of 
(smooth) metrics on Cq. While one can always consider constant curvature, here flat, 
metrics with geodesic boundary curves, it turns out that this particular representative is 
in general not the natural one to flow surfaces with boundary towards minimal surfaces. 
In particular, one would like to avoid the possibility that a boundary curve of the domain 
(on which we after all impose our boundary condition) can shrink to a point and thus 
be lost. 

For the cylinder we shall thus consider smooth metrics compatible with c which have 
constant curvature —1 and for which the boundary curves have both the same constant 
geodesic curvature. 

We first recall the following standard fact of complex analysis 

Lemma 2.1. To any smooth conformal structure c on Co there exists a unique number 
Y > 0 such that (Co, c) is conformally equivalent to ([—Y,Y] x 5' 1 ,ds 2 + d9 2 ). 
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On such a cylinder ([—Y, Y] x S 1 ,ds 2 + d6 2 ) we can then use the following hyperbolic 
metrics whose structure is well known from the Collar lemma [11] 

Lemma 2.2. On ([—Y, Y] x S 1 ,ds 2 + d9 2 ) there is a one parameter family of collar 
metrics 

gi = pe{s) 2 (ds 2 + d8 2 ) 

where 

Pi(s) = o-“7777’ £ e (0, L 0 (Y)), L 0 := 

27TC0S(^;) Y 

which are all hyperbolic and whose boundary curves have the same constant geodesic 
curvature k = nty ■ 


As admissible metrics for our flow we shall thus consider 

A4_i := {f*g : / : Co —> [—Y. Y] x S 1 smooth diffeomorphism, 

g = ge a collar metric as in Lemma 2.2, £, Y £ (0,oo)}. 


( 2 . 1 ) 


While Lemma 2.2 does not yet give a canonical representative of a conformal class, such 
a choice can be made so that the following splitting of the tangent space is respected 

Lemma 2.3. For any g £ A4_ 1 we have 

T g M _i = {L x g : X £ r(rC 0 )} © Re{U{g)) 0 span{^ ■ g}. 

where Re(fH(g)} is L 2 (Co, g)-orthogonal to {Lxg} © span{if 2 ■ g}. 

Here H(g) is the real vector space of quadratic differentials that are holomorphic in the 
interior of Co, continuous upto the boundary and whose traces on dCo are real. Further¬ 
more r(TCo) stands for the space of smooth vectorfields on Co which are tangential to 
dCo on dCo and if g :Co^-S.is characterised by 

*p 2 g -9 = r(i(p 2 e )(ds 2 + d0 2 )) 

for g = f*ge 0 , gi the collar metrics of Lemma 2.2. 

We recall that for the cylinder the space H is simply made up by elements of the form 
cdz 2 , c € R, z = s + i9, for collar coordinates (s, 6) £ [—Y, Y] x S 1 as in Lemma 2.2. 

We also remark that the orthogonality relation claimed in the lemma is a simple conse¬ 
quence of the fact that the real part of a holomorphic quadratic differential is trace and 
divergence free. 

This lemma implies that the most efficient way (i.e. with least L 2 velocity) to lift a curve 
[<?(•)] from Teichmiiller space A4~i/T>o to M_ 1 is as a horizontal curve, moving only in 
the direction of Re('H(g)). Here T>o denotes the space of smooth diffeomorphisms from 
Cq to itself that are homotopic to the identity. 

For cylinders we can describe such horizontal curves of metrics explicitly by the following 
lemma which is proved in the appendix 

Lemma 2.4. Let g > 0 be any fixed number. We define Y = Y V : (0, oo) — > (0, oo) by 

Y(f) = j(| — atanfq ■ £)) 
and f e = ft : C 0 [-¥{£), Y(£)\ x S 1 by 

fl{x,<t>) = («? 0 * 0 . 0 ) = {jfatan(^- ■ tan(|£a;)), f>) 
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where £o = * s determined through the condition Y(£q) = 1- 

Then the family Gz := (fe)*[p'j{s)(ds 2 + dO 2 )) is horizontal, i.e. for every t 

^G t eRe(H{C 0 ,G t )). 

Since Y(-) is a bijection, we can combine the above result with Lemmas 2.2 and 2.1 to 
conclude that any horizontal curve of metrics in A4_i must be of the form /*(G^) for 
some fixed ij > 0 and a fixed diffeomorphism / : Cq — » Cq. 

As such, we shall from now on consider rj > 0 to be fixed and will in particular allow 
all constants to depend on this number as well as on the boundary curves T± (and their 
parametrisations a±) without further mentioning this. 

To describe the space of admissible maps, we first recall that the prescribed boundary 
curves T± are assumed to be disjoint, closed C 3 Jordan curves of which we shall fix 
proper C 3 parametrisations 

a± : S 1 —> T±. 


We then consider maps in the space 

H^(Co,g) '■= {u £ LJ 1 ((Co, 5 ),K") such that u : dC± —> r± is weakly monotone } 
i.e. H 1 maps so that the traces u\ac± can be written in the form 

u\dC ± = a± o ip± 

for some weakly monotone functions ip± : S 1 —> S 1 . Here and in the following we identify 
dC± := {±1} x S 1 with S 1 when convenient. 

It is well known that the space is not closed under weak H 1 convergence as one 
can find sequences of maps with bounded energy for which the boundary curves T± are 
parametrised over smaller and smaller arcs of dC± thus resulting in a weak limit that no 
longer spans T±. The standard way to deal with this loss of completeness is to impose a 
three-point-condition. So we shall restrict the set of admissible maps for our flow to 

Hr *(Co) := {« £ ^r(Co) : u\ dc± = a± o <p± for <p± satisfying 

2 tt ( 2 . 2 ) 

<P±(8k) = 9k for dk = —k, k = 0,1, 2}. 

To compensate for the (in our case 6 ) lost degrees of freedom we need to allow the 
metric to move not only in horizontal direction but also through the pull-back by select 
diffeomorphisms. 

For this purpose we will define (and discuss) a suitable family of diffeomorphism hb^, 
<j> = (</>+, (f~) £ R 2 , b = (b + ,b~) £ C 2 , with l^l < 1 later on in section 4.1.1. We remark 
that by using the one Killing field that is available for the cylinder we could reduce the 
number of degrees of freedom to 5 instead of 6 (e.g. by asking that (f> + + fr = 0) though 
this would not lead to a significant simplification. 

All in all we then say that a curve ( u,g)(t ) is an admissible representative of a curve of 
equivalence classes [(u,g)(t)] £ A if 

• g{t) = ^ 0 (t) b(t)dW f° r a horizontal curve of metrics g(-) £ A4_i and continuous 
families of parameters (b, (/))(■) £ where flh := (l?i(0 )) 2 x I 2 C C 2 x l 2 is the 
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domain of parameters for the diffeomorphisms hb t< j> that will be defined in section 

4.1.1. 

• u(t) £ H^ ^Cq, R n ) for every t. 


We remark that we can and will assume without loss of generality that the initial metric 
g(0) of the horizontal curve ( g ) is given by one of the metrics Gi described in Lemma 
2.4 simply by pulling-back the whole problem (including the parametrisations a± used 
in the three-point-condition) by a fixed diffeomorphism. As such we shall from now on 
consider metric in the set 

M. := {hl ^Ge. : ( b,4 >) £ fi/j, i £ (0, cxo)}. 


2.2 Definition of the flow 

As we consider a problem with a Plateau boundary condition, the space of admissible 
variations does not form a vectorspace. As such the flow that we shall define will not be 
governed by a system of PDEs with prescribed boundary values but rather, for the map 
component, by a partial differential inequality. 

To motivate the following definition we first make some general formal computations, 
which we of course do not claim to be new in any way, but which are rather included for 
the convenience of the reader. We remark in particular that the differential inequalities 
we derive correspond to the ones obtained in [1] and [6] in case of the domain being a 
disc. 

Given a functional T defined on some (Hilbert)manifold B we may want to define a 
gradient flow under the restriction that the velocity dtw at each time is constrained to 
some closed convex cone X(w(t)) C T w mB, e.g. because we want to constrain the flow 
to some convex set A and thus the velocity to the corresponding solid tangent cone. 

One can formally define such a gradient flow by asking that 

d t w = P x{w) ( - VJ(w)) (2.3) 


where P X W : T p B —> X(p) is the nearest point projection. 

We then observe that a variational formulation can be given by asking that at each time 
the velocity dtw is given by a variation -^\ E — 0 w e of wo = w(t) which is admissible in that 
-j;\e=oUJ E £ X(w(t)) and which, among all such variations, minimises the functional 

f e \ E=0 T(w E ) + ^\\£\ EM \\ 2 . (2.4) 


In practice, such a formulation often asks for more regularity, in particular of dtw , than 
what we can a priori expect of a weak solution. So consider instead 2-parameter families 
w E} s with w E: o = w E such that each of the families 1 u. t s gives again an admissible variation 
of w(t). Then if w E p = w E minimises (2.4) we must have that 

Zs\s=o-^\s=oX(w Et s) + (-jg\s—o-^\ E= ow E ^,dtw) > 0 (2-5) 

which gives not only a weaker condition than (2.4) but often requires less regularity of 
dtw than (2.4). 
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Going back to our problem of defining a gradient flow of the Dirichlet energy on the 
set A we recall that the (negative) L 2 -gradient of the energy with respect to the metric 
variable can be written in the form jRe(<&(u : g)), where &(u,g) is the Hopf-differential 
which is given in isothermal coordinates (s, 0) = {s g: 0 g ) of ( Co,g ) by $(u,g) = (|u s | 2 — 
\ug\ 2 - 2i{u s , ug)) ■ dz 2 , z = z g = (s + id). 

The weak formulation (2.5) thus translates to the condition that with v := -^g\s=e=oU E ,s 
and h = -j^[s\s=s=oge,6 

{dv, du) g + v ■ dtu dv g 

for all variations {u e ^,g E .s) of map and metric that are admissible in the sense described 
above. 

On the one hand, the resulting differential inequality for g can be simply recast as a 
differential equation 

d t g=-^{Re{^{u,g)) (2.7) 

to be solved on AT Here is the L 2 -orthogonal projection onto the tangent space V(<?) 
of M = {hi : {b, tf>) £ f l £ ( 0 , oo)} which is given by 



(--Re($(u,flO) 


dtg , h) dv g > 0 


( 2 . 6 ) 


V{g) := Re{H{g)) © {L h ^ x g : X £ X(b, </>)}, (2.8) 

for g = hl ^Gi where X(b,<j>) is the 6 dimensional space of vectorhelds generating the 
diffeomorphisms kg compare section 4.1.1. 

On the other hand, admissible variations of the map can be described as follows. Given 
u £ ffp *(Cb) we let tp± be such that u\ac± = ck± ° <p±- Then functions of the form 
a± o (<p± + £ ■ (3± + 0(e 2 )) are again monotone parametrisations of T± at least for e in a 
small onesided interval [0,£o) if /3± can be written in the form j3± = A± ■ (ip± — ip±) for 
some numbers A± > 0 and weakly monotone functions ip±. 

As variations ^| E= o u of the map component we thus consider elements of 

T u H t,*( c o) ■= {v G fl , 1 (C' 0 ,M") :v\ dc± = A ±a'±{<p) • (ip± - ip±) for A± > 0 and 

V>± G C°{S 1 ,S 1 ) weakly monotone with ^±{0k) = Ok}- 

(2.9) 

We remark in particular that to any v G T+iJp *(Co) there is a onesided variation 
(u E ) C Hp ^(Co), £ G [0,e- 0 ) with -£\ e=0 u e = v, see [ 6 , Lemma 2.1]. 

As r+m „ is in general not a vectorspace, but only a convex cone, we cannot reduce the 
resulting partial differential inequality 

{du,dw) L 2(C 0 ,g) + / w ■ dtudv g > 0 for all w G T+ *(Cb) (2-10) 

to a PDE with a standard boundary condition though one immediately obtains that u 
satisfies the heat equation d t u = A g u in the interior. 

Furthermore, as pointed out in [ 6 ], the additional condition that 

/ {du, dw) + A gU ■ w dvg > 0 for all w G T+iLp *, 
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can be seen as a weak Neumann-type boundary condition. 

Given that (2.7) and (2.10) were motivated by the idea that dt(u,g) should minimise the 
functional (2.4), the so called stationarity condition, asking that 


1 

4 


Re($(u,g))L x gdv g 


Du(X) ■ A gudvg = 0 for every X £ T(TCq )* (2-11) 


where 

r(TC 0 )* := {Y £ r(TC 0 ) : Y(± l,0 fc ) = 0}, 

results if one considers variations of the form [u(t + e) o f e ,g). 

Similarly one expects the energy to be non-increasing along the flow, compare (2.14) 
below. 


All in all, we define 

Definition 2.5. A weak solution of Teichmuller harmonic map flow on the cylinder Co 
is represented by a curve of maps 


u £ L°°([0, T), H^{Co, R")) n ^([0, T) x Co) 
and a curve of metrics g £ C 0,1 ([0, T), M) which satisfy 

/ (du, dw) g + d t u ■ w dv g (t) dt > 0 for all w £ L 2 {{ 0, T], Tff *(Co)) 

J[0,T]xC o 

and 

dtg = \Pg {Re($(u,g)) for a.e. t. 

Such a weak solution is called stationary if it satisfies (2.11) for almost every t, 
say (it, g) satisfies the energy inequality if for almost every t\ < t 2 

E{u,g){ti)-E{u,g){t 2 )>\[ [ \d t u\ 2 dv g dt+ f \\d t g\\ 2 L 2 {c ^ g) dt. 

Jti J Co " 


( 2 . 12 ) 

(2.13) 
and we 


(2.14) 


2.3 Main results 

For the flow we just defined we will prove the following two main results 

Theorem 2.6 (Existence of global solutions). Let T± be two disjoint closed C 3 Jordan 
curves. Then to any initial data (uo,go) £ *(Co) x A4 there exists a stationary weak 

solution (it, g) of Teichmuller harmonic map flow which is defined for all times, smooth 
in the interior of Co and satisfies the energy inequality. 

The above solution flows to a minimal surface in the sense that 

Theorem 2.7 (Asymptotics). Let ( u,g){t ), t £ [0,oo), be a stationary weak solution 
of Teichmuller harmonic map flow that satisfies the energy inequality and for which the 
three-point-condition does not degenerate in the sense that limsupj^^ 1 — |6 ± (t)| > 0. 
Then there is a sequence of times ti —> oo such that the equivalence classes [(it, <?)(tj)] 
converge to a critical point of the area in one of the following ways: 


I (Non-degenerate case) If inj(Co, g(U)) 0 for i —> oo then f*(u(ti), gitf)) con¬ 

verges to a limit (moo,(?oo) where g ^ G Ai and where Uoo G *(Co,K™) DC 0 (Co) 
is a (possibly branched) minimal immersion. 

Here 

fi ■■= ho, 2 n ni , nf= L^J (2.15) 

and the convergence for the metric component is smooth convergence on all of Cq 
while the maps converge uniformly on the whole cylinder Cq as well as strongly in 
i? 1 (Co) and weakly in H( oc (Cq \ (Jj±Pj Z ) away from the points P± = (±1 ,0j), 
j = 0 , 1,2 at which the three-point-condition is imposed. 

II (Degenerate case) If inj(Co, g(U)) —> 0 then f* (u(U), g(ti)) converges locally on 
Co \ ({0} x S 1 ) to a limit (uoo,Soo) which is such that 

• each of the cylinders (C±, c/oo), C± := {0 < ±s < 1} x S 1 is isometric to the 
hyperbolic cusp 

([ 0 , oo) x S 1 , p 0 (s) 2 (ds 2 + d0 2 )), p 0 (s) = 1 

2ng + s 

• The two maps Moo|c± can be extended across the punctures to give two (possibly 
branched) minimal immersion u ± G H( ± *(D) D C°(D) parametrised over 
closed disc in R 2 each of which spans the corresponding boundary curve . 

Here the convergence is smooth local convergence for the metrics and weak Hf oc 
convergence on (C_ U C+) \ (J ■ l ±Pf L as well as locally uniform and strong H) oc 
convergence on C_ U C+ for the maps and the diffeomorphisms are again given by 
(2.15). 


We remark that while we only obtain convergence in H 1 fl C° respectively in H 2 away 
from Pp, the limit u 00 is indeed far more regular than that. Namely, classical regularity 
theory for solutions of the Plateau-Problem, see e.g. [17] or [5], yields that Uoo is of class 
C 2 ’“, a < 1, upto the boundary. 


3 Short-time existence of solutions 


We shall prove short-time existence of solutions to arbitrary initial data based on a time 
discretisation scheme. We remark that this method has been carried out successfully to 
obtain solutions of several other geometric flows, e.g. by Haga et. al. [ 8 ] for harmonic 
map flow and by Moser [10] for biharmonic map flow, and that also the solutions for 
the evolution to minimal discs by Chang-Liu [1] respectively to discs of prescribed mean 
curvature of Duzaar-Scheven [ 6 ] were obtained this way. 

A key part of this section consists in proving suitable a priori estimates for the approxi¬ 
mate solutions resulting from such a time discretisation. For some of these estimates we 
will be able to appeal to work of Duzaar and Scheven [ 6 ] whose delicate estimates allowed 
them to prove H 2 bounds upto the boundary but away from the points Pp despite their 
equation being non-linear. 

In the present paper the challenges are somewhat different as we do not have to deal 
with a non-linear equation for the map component but instead have to understand the 
interplay of the map and the metric component of the flow. What makes this particular 
aspect of the flow quite delicate, is that this relation involves a non-local projection 
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operator. This forces us to prove estimates that are valid not just near most boundary 
points but rather in neighbourhood of every boundary point, including the points Pp at 
which we impose the three-point-condition. 


3.1 The time discretisation scheme 

To begin with we outline the time-discretisation scheme and show that it is well defined. 

Given an initial pair (ito, 5 o) £ Hr ♦(C’o) x At and a (small) number h > 0 we construct 
an approximate solution of Teichmiiller harmonic map flow using the following time- 
discretisation: 

For j = 0,1, 2,.. we let tj = t'f = j ■ h, set u h (t) = uq for t £ [to, ti] and then construct 
iteratively the approximate solution ( u h ,g h ) on the interval [tj,tj +1 ] as follows: 

First determine g h {-) on (tj,t^ +1 ] as the solution of 

dtg(t) := ^pv {t) ($(u h (tj),g(t)) with g{tj) = g h (tj). (3.1) 

Then select u h (tj+{) as a minimiser of the functional Pg h ^ t . +1 ^ u h where 

Hg, v (w) = E{w,g) + ^r\\w - v\\ 2 LHCotg) . (3.2) 


The existence of a minimiser of this functional is assured by the direct method of cal¬ 
culus of variation thanks to the H 1 - weak-lower semicontinuity of u K > E(u, g) as well as 
the Courant Lebesgue Lemma and the resulting equicontinuity of the traces u\qc 0 i c.f. 
appendix A.l. 


To be more precise, we have 

Lemma 3.1. For any g £ A4_i, any map u £ F[f t (Cb) ani ^ an y ^ > 0 there exists a 
minimiser w £ Flf ^(Co) of 

T(w) =E{w,g) + jk\\w - u\\ 2 L 2 (Cotg) 

and w satisfies 



u ) • vdvg > 0 for all v £ Tff *(Co) 


in particular A g w = ^ • (w — u) in the interior of Cq. 
Furthermore w satisfies the stationarity equation 


1 

4 


J Re($(w,g)) ■ Lxgdvg + J dw(X) ■ A g w dv g = 


0 for all X £ r(TC 0 )* 


and the energy inequality 


E (w , g) + jh'\\w- u\\l^ Co , g ) < E(u , g). 


(3.3) 


(3.4) 


(3.5) 
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We remark furthermore that the minimiser w is bounded by 

IMU“(C 0 ) < ||«||i-(Co) ( 3 - 6 ) 

so that the L°° norm of the map component of the flow is non-increasing in time. Indeed, 
if the above estimate would not be satisfied, we could compose w with the nearest point 
projection to the ball with radius ||it||z,=° to obtain a function with smaller energy J-. 

In the present setting of metrics on a cylinder, short-time existence of a solution to 
the differential equation (3.1) on M. is a simple consequence of the fact that ji is one¬ 
dimensional since this means that the evolution of the metric could be expressed as a 
system of (in total 7) ordinary differential equations. As such it is easy to check that the 
projection satisfies the following Lipschitz- estimates 

Lemma 3.2. Let K be a compact subset of the set of admissible metrics A 4 = {g = 
£ (0,oo), (b,(f>) £ flh} and let be the L 2 -orthogonal projection onto V(g) := 

T g A4. Then 

||pV(i2 e ($ 1 ))_pv ( ^( $a ))|| cfc < C || 5l _^|| c ,.||$i|| L1 ( Co)+ q| ite (^ 1 )- jRe ($ 2 )|| L1(Co) 

for all gi £ K, quadratic differentials 4q on (Co, <?i), i = 1, 2 and k £ N, where C depends 
only on k and K. Here and in the following the C k norms are computed with respect to 
a fixed coordinate chart. 


We remark that the real part of the Hopf-differential can be written equivalently as 

Re($(u, g)) = 2 u*g^ - \du\ 2 g ■ g (3.7) 

so that we can bound the differences of Hopf-differentials by 

||i?e(d>(u, 5 )) - Re{$(u,g))\\ L i < C\\g - g\\ C o +C\\u-u\\ H i (3.8) 

for a constant C that depends only on a bound for the energies E(u 1 g) 1 E(u,g). 

Combined, Lemmas 3.1 and 3.2 thus imply that solutions of the time discretisation 
scheme exist for as long as the injectivity radius inj(Co,g) = 2H of the domain (Co, g) is 
bounded away from zero and infinity and the parameters (6, <p) remain in a compact set 

of fl h - 

We furthermore remark that the energy of such an approximate solution is non-increasing, 
namely on the open interval (tk 7 tk+ i) it decreases by 

E{{u h , g h ){t k )) - E((u h (t k ),g h (t k+ 1 )) = f \\dtg h \\ 2 L 2 { c 0t gh)dt 

J [tkJk+ 1 ] 

= 4 f \\ P v ( Remg h ( t ), u h ( t k Mh iCo ^)dt 

10 J[t k ,tk+l] 

while at t k+ 1 there is a further loss of energy of no less than 

E(u h (t k ),g h {t k+ 1 )) - E{{u h 7 g h ){t k+1 )) > ^-|| u h {t k ) - u h {t k+ i)\\ 2 L 2 (Cotgh{t K +i)) , 
compare (3.5). 

All in all we can thus estimate 

E((u h ,g h )(t k )) < E((u\g h )(t- k )) - f" 111^11^,^) + 11^111^,^ ( 3 -9) 

J tz. 
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for k < k where compute the norm of the difference quotient 

D b u h (x, t) = jj • ( u h (x , t + h) — u h (x, t)) 
with respect to the piecewise constant curve of metrics g h \[t k ,t k+1 ) = 9 h (tk+i)- 
We obtain in particular that the length of the curve of metrics is bounded uniformly by 
L L 2 (g h \[o,t]) < E^ /2 t 1/2 for every h > 0, 

Eq an upper bound on the initial energy E(u(0),g(0)). 

Consequently, to any given (uo,go) there exist numbers e > 0, T > 0 and C < oo such 
that for every parameter h > 0 the solution of the time-discretisation scheme exists 
at least on the interval [0, T] and so that the metric component g h = h b ^Gt satisfies 
estimates of the form 

inj(Co, = 2£> e and \b±\ < 1 — e, |^ ± | < C (3.10) 

on this interval. 

Similarly we could obtain an upper bound for £ on [0, T], but we remark that £ is indeed 
bounded from above uniformly in time in terms of only the initial energy. Namely, 
let dp := dist(r + ,r_) > 0 be the distance between our prescribed disjoint boundary 
curves. Then the energy of any map w £ H^(Cq) with respect to a metric g of the form 
g = 0 Gn is bounded from below by 

i r 2 " r YW 2 2 

E(w,h* b M e ) = E(w,G t ) = - / Kl 2 + \wo\ 2 dsde 

1 Jo J-Y(e ) 



for some c > 0 and w = w o h b ^. The resulting lower bound on Y (£) results in an upper 
bound for £ and thus the injectivity radius of the form 

£ < L(E 0 ), (3.11) 

where L(Eq) depends only on an upper bound Eq on the initial energy E(uQ,go) (and as 
usual the geometric setting and the fixed number rf). 

This uniform upper bound on £ and the resulting control on the metrics Gt near the 
boundary 8Cq will allow us to derive a priori bounds for the map component near 8Cq 
that are independent of inj(Co, g). This will be crucial for the asymptotic analysis (in the 
degenerate case) carried out later on. Conversely, all estimates near the boundary will 
depend on b as the case | 6 ±| —> 1 corresponds to the degeneration (in collar coordinates) 
of the three-point-condition. 

To prove that the above time-discretisation scheme converges to a solution of the flow, 
we need to derive a priori estimates for the map component where we will distinguish 
between 

• the interior of the cylinder where standard estimates for the heat equation apply 

• the boundary region away from the points = (±1 ,0k), on which we shall be 
able to appeal to results of Duzaar-Scheven [ 6 ] 

• the region near the points at which the three-point-condition is imposed. 
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3.2 A priori estimates in the interior and near general boundary 
points 

Let u be any fixed map, h > 0, let g = h J ^Gi be a metric as in Lemma 2.4 and let u be 
a minimiser of the functional F gU = E(u, g) + j^\\u — u\\ 2 L2 ^ Co g )- 

Then, setting f = L(u — u) we know that 

f ( du,dv)dv g + f v ■ fdvg > 0 for v £ T^H^ .*(Co)- (3-12) 

J Co J Co 

In particular A g u = f in the interior, so standard elliptic estimates combined with the 
upper bound (3.11) on i yield 

Lemma 3.3. Given any to > 0, Eo < oo and any 8 > 0 there exists a constant C < oo 
such that the following holds true. Let u £ „(Cb) be any map of energy E(u,g) < Eq 

which satisfies (3.12) for some f £ L 2 (Cq) and g £ At with 2inj(C'o,<?) > £q. Then 

|V >| 2 + \duf g dv g < C • (E(u,g) + ||/||| 2(CoiS) ). (3.13) 

We remark that the region in which Ge degenerates as £ —> 0 is contained in what 
corresponds to arbitrarily small cylinders [—<5, <5] x S 1 with respect to the fixed coordinates 
(x, 6 ) £ Cq = [—1,1] x S 1 since we use hyperbolic rather than flat metrics to represent a 
conformal class. Therefore 

Remark 3.4. The analogue of (3.13) is valid with a constant independent of on every 
compact region of (—1, 1 ) x S 1 \ ({0} x S 1 ) and for every metric g £ A4 as well as for 
metrics hb^Gi— o described in (A. 13). 

Near the boundary but away from the points Pp we can use the results of Duzaar- 
Scheven [ 6 ] which apply to more general (in particular non-linear) equations. Namely, as 
explained in appendix A.l, we can derive the following a priori estimates from Theorem 

8.3 of [ 6 ] 

Proposition 3.5. For any bo < 1, £q > 0 and Eq < oo there exist constants £i,rg > 0 
and C < oo such that the following holds true. Let g = ^Ge for some l> to and 
|^ )± | < bo. Suppose furthermore that f £ L 2 {Co,g) and that u £ i/p t (Co) has energy 
E(u,g) < Eq. Then, if u is so that (3.12) is satisfied for all variations v £ T + iLp *(Co) 
with support in a ball Bf.(jp), where p £ Co and r £ (0,ro) are such that 

B 3 r (p) n {Pf} = 0 

and if the energy on this ball is small in the sense that 

E(u,B 3 (p)) := \ f \du\ 2 dv g < £\ 

Jb°( p ) 

then u £ H 2 (B 3 ^ 2 (p),g) with 

[ \V 2 g u\ 2 + \\7 g ufdv g < ^E(u,B 3 (p))+C [ \f\ 2 dv g . (3.14) 

J B a r/2 {p) r Jb°( p ) 


'[—1+5,1—5] xS 1 
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Here and in the following we denote geodesic balls in (Co, g) by B g (p) := {p £ Co : 
d g (p, p) < r} and compute the energy on a geodesic ball with respect to the corresponding 
metric unless indicated otherwise. 

As we shall use this and the subsequent lemmas to control the map near the boundary, 
it is important to remark 

Remark 3.6. The above result is valid also without imposing a lower bound on £, and 
in particular also for the metric Go defined in (A. 13), as long as one considers only points 
contained in a compact subset K C Co\{0} x S 1 and allows the constants to depend also 
on this set K. Similarly, on compact sets K CC Co \ ({0} x 5 1 ) the a priori estimates 
derived in the subsequent Lemma 3.7 and in Corollary 3.12 are all valid for metrics h £ ^Gt 
with l > 0 and | 6 ±| < bg < 1, again with a constant that also depends on K. 


As our target is Euclidean space which ’supports no bubbles’, i.e. for which there are no 
non-trivial harmonic maps from S' 2 , we can furthermore exclude a concentration of the 
energy near general points of the boundary 

Lemma 3.7. To any numbers A, Eg < oo, bg < 1, d, tg > 0 and any e > 0 there exists 
a radius r > 0 such that the following holds true. 


Let u £ Hf *(C 0 ) be a map of energy E(u,g) < Eg which satisfies (3.12) for a function 
f £ L 2 (Co,<?) with II/Hl 2 < A, a metric g = h^Ge, with £ > £g and | 6 ± | < bg and 
variations v £ T+iJp *(Cq) with supp{v) C C* := Cg \ (J y ± . 


Then the energy is small 
on balls around arbitrary points 


E{u,B g r (p )) < e 


(3.15) 


p £ C*(d) = C* g (d) := Co \ 1J B B d (lf). 

i>± 


In particular, the estimate 

IMI H 2 {c*(d)) < C • E(u, g) + C ||/|| 2 2 
holds true with a constant C that depends only on A,Eg,£g, bg and d. 


Proof. In order to prove the first part of the lemma we argue by contradiction. So assume 
that for some numbers e,d > 0, foo < 1 and Eg, A < oo there is a sequence of (■ Ut,gi,fi ) 
as in the lemma and a sequence of radii f, — > 0 for which sup^, ec;.(d) E ( u > > e - 

Here we can of course assume that e < £i, the number of Proposition 3.5. 

We first prove 

Claim: There exist radii Vi —> 0, points pi £ C*. (d/2) and numbers Xi —> oo so that 

E(u, B gi (pi)) = s = max E(u,B gi (p)). 

P eB 


To prove this claim let us first choose points pi and radii vy —> 0 so that E(u, B gi (yi)) = 
e = maxpgc 1 . ( d ) E{u, B-jf (p)). Then the claim is trivially true for pi = y, unless the points 
Hi converge to the boundary (relative to Co) of the set C*.{d) defined in the lemma. 
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So assume that, after passing to a subsequence, dist 9i (y l . P ±) —> d for one of the point 
Pf, say for P+. 

We then consider concentric annuli 

A i'-= Bf_ 2kri \B^_ 2{k+1)ri (P+) 

constructed so that two balls of radius one having its centre in A & the other in A /.+2 
are always disjoint. Thus the number of such annuli that contain a point p for which 
E(u , B r f ( p )) > e can be no more than K e = 2 \_Tf-\ . 

We now choose Ni —> 00 so that for i large K e ■ (W; + l)ri < d/2 and observe that 
B\ B^ 2 (x 0) must contain an annulus A i °f thickness N^n which does not 

contain any point p with E(u, B/* (j>)) > e. Possibly reducing r t so that the maximum 
of p 1 —y E(u, B-/\ (p)) on B 2d \ B d _ k r . (xo) is equal to e and selecting pi to be a point at 
which this maximum is achieved then implies the claim. 

Based on this claim we now derive a contradiction using a standard blow-up argument 
where we distinguish between 

Case 1: dist fli (pj, dCo)r~ l —> 00 for some subsequence 
Case 2: dist(p i; dC^r// 1 < C for some constant C < 00 

In the first case, we rescale the maps to maps Vi(x) = u(exp//. (r^x)) that are defined on 
larger and larger balls in R 2 . We also observe that as always in such a bubbling argument 
the resulting metrics can be written as (exp®* (ry))*<7i = rfgt for metrics gi that converge 
locally to the euclidean metric g euc i- 

In particular, the H 2 bounds on subsets of (Co, gi) obtained in the previous lemmas imply 
that the iP-norm of the maps u,; on compact subsets of (R 2 , g euc i ) are bounded uniformly 
and that A Vi —>• 0 locally in L 2 since we have assumed ||A ffi Ui|| i 2 to be bounded. 

Thus, after passing to a subsequence, we conclude that Vi converges strongly in H/ oc and 
weakly in H/ oc to a limit Vao : R 2 —> R n which is harmonic and has bounded energy 
and is thus constant. At the same time P(uoo,Pi(0)) = lim^oo E(vi, -Di(O), g/) = 
lim^oo E(ui, B/*(pi)) > e resulting in the desired contradiction. 

In the second case we rescale not around the points pj themselves, but rather around 
their nearest point projection pj to the boundary of Co- The resulting maps, defined 
on larger and larger subsets of the halfplane H = {y £ R 2 , j/ 2 > 0}, satisfy uniform H 2 
bounds on compact sets of H and have energy at least e in the ball Bc+ 1(0) fl H. 

We again obtain a harmonic limit with bounded energy, now defined only on the half¬ 
plane, but furthermore constant on the axis cffl, since the maps Ui\gc 0 are equicontinous, 
compare Corollary A.2 in the appendix. Thus also this limit must be constant leading 
again to a contradiction. This completes the proof of the first part of the lemma. 

The second part now immediately follows from Proposition 3.5 and the first part if we 
choose e = e 1 to be the constant of that proposition. □ 

A further consequence of Proposition 3.5 is 

Corollary 3.8. Let u € H p *(Co) be a function for which (3.12) is satisfied for some 
f € L 2 (Co) and g = hi £ > 0, and for all variations v € T+ H/ * (Co) with supp{v) C 
C* := Co \ U; ± Pj 1 ■ Then u is continuous on all of Co, in particular in the points 
where the three-point-condition is imposed. 


15 


Furthermore, let u,i £ Hf t (Co) be any sequence of maps with uniformly bounded energy 
for which (3.12) is satisfied for functions fi with sup, \\fi\\L 2 (M,gi) < oo and metrics gi = 
hl. ^G lt with sup 4 16^1 < 1 and infj^j > 0 and again for variations v £ T+iJp *(Co) 
with suppiv ) cC*. If no energy concentrates at the points P± in the sense that 

lim sup E(ui, B 9i (Pf 2 )) = 0 j = 0,1,2 (3.16) 

r —>0 j J 

then the maps Ui are equicontinuous on all of Co. 

Proof. Let Ui be a sequence of maps as described in the lemma. As Lemma 3.7 yields 
uniform Lf 2 -bounds and thus equicontinuity for the uf s on any compact subset of C* it 
is enough to prove that to any number e > 0 there exists a radius ro so that 

OSC Ui < e, j = 0,1, 2. 

B g r' 0 (P±) 


To begin with, we recall from Corollary A.2 that the traces rti|ac 0 are equicontinuous so 
that for ro sufficiently small 

osc Ui < sup osc Ui + osc Ui < sup osc Ui + he (3-17) 

re(0,r o \dBB^ndGo re(0,r o ] 


where dB r C Co is the boundary relative to Co and where we consider balls with centre 
Pj unless indicated otherwise. 


We bound the oscillation over dB r by deriving suitable iL 2 -estimates on annuli A* := 
B 9 s \ B I' . First of all, by (3.16) we can assume that ro is small enough so that 

4 r 4 r 


E(ui, B^o) < £1, the number of Proposition 3.5. Then, given any r G (0,ro] we cover 
the above annulus A l r by balls B 9 C(xk) so that the corresponding balls with double 

the radius are contained in \ {Pj^} and so that no point is contained in more than 
K of these larger balls B 9 \ 2 (xk), K independent of r and i. Since (3.3) is satisfied for 
(ui, gi, fi), at least for variations supported on C *, we can apply Proposition 3.5 to bound 


r / 4 ( x fe) 

and thus also 


[ \V 2 gi Ui \ 2 dv gi <Cr- 2 E( Ul ,B a r Ux k ))+C\\M \ 2 Lhb h , Xk)) (3.18) 

JB g j.Ax k ) V k ” 


J |V^i| 2 dv gt < Cr~ 2 E(ui, P|;) + C\\M 2 L2(Bt) (3.19) 


for a constant C that is independent of i. 


Observe that while the oscillation is invariant under the rescaling fq( x) = Ui(exp p ± ( rx )), 

J 

the left-hand-side of the above estimate transforms as 


/ |VL>| 2 ^ < Cr 2 f | W 2 giUi f dv gi + CE(ui, Pf*) 

J A x ( 0) J A r 

< CE(ui, B%). o ) + Cr 2 ||/i|| 2 2(B 9^ } . 


(3.20) 


Applying the Sobolev embedding theorem on a suitable subset of the fixed half-annulus 
Ai(0) C R 2 thus allows us to conclude that 


( osc uA 2 < C ■ 

v dB gi 



+ \dui\ 2 dx < CE(ui,B^ o ) + Cr 2 \\fi\\ 2 L 2 i B^ o ) 
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for every r £ (0,ro] again with constants that are independent of i. 

Since the L 2 -norms of fi are uniformly bounded and since we assumed that there is no 
concentration of energy at the points P± we can thus choose ro small enough so that the 
above expression is less than e/2 as desired. □ 

Remark 3.9. We observe that the claim of the above corollary remains true on arbitrary 
compact regions of Co \ ({0} x S 1 ) also without the assumption of a lower bound on t. 
Similarly, as all arguments are carried out locally, knowing that u satisfies (3.3) for 
variations supported in a set U\{ J ■ ± Pj is sufficient to conclude that u is continous on 
every compact subset K of U with modulus of continuity depending only on K and the 
bounds on |6 ± |, E(u,g) as well as the local L 2 -norm of /. 

3.3 No concentration of energy at points Pp. 

As the flow of metrics is determined in terms of the (non-local) projection of the Hopf- 
differential onto V, we need to exclude the possibility that a non-trivial amount of energy 
(and thus possibly of L 1 norm of <f>) is concentrating near one of the points Pp. Such 
a concentration of energy would be lost in a limiting process meaning that we could not 
expect that the evolution of the limiting metric would be described by the projection of 
the limiting Hopf-differential. 

To this end we prove the following key lemma 

Lemma 3.10. To any given numbers A,M,Eo < oo, bo < 1 and £,do > 0 there exists a 
radius r > 0 such that the following holds true. Let u £ H / *(Co) be a map with energy 
E(u,g) < Eq that is bounded by ||u||l<*> < M and that weakly solves the differential 
inequality (3.12) for a metric g = hf ^Ge for which < bo and for a function f with 
ll/llz, 2 (Co,s) — A- Then the estimate 

E(u,B 9 (x 0 )) < e 

holds true for every point Xo = (x, 8) £ Co with |x| > do, in particular for Xo = Pf. 

Proof. Thanks to Lemma 3.7 it is sufficient to establish the claim for the points P- , say 
for xo = P 0 + . 

So assume that the claim is wrong because for some fixed numbers A, M, Eo < oo and 
bo < 1 there exists a number £2 > 0 such that there are triples (it*, gi = h'f ^.G^, fi) for 
which all the assumptions of the lemma are satisfied, but for which energy concentrates 
at xo in the sense that 

E(ui,B 9 i(x 0 )) > £ 2 

for a sequence of radii r* —> 0. 

We remark that the cliffeomorphisms hb^ defined later on in section 4.1.1 are such that 
hb^ = h b ^ in a neighbourhood of dCo if the parameters and dA agree modulo 27 t. 
Thus, after passing to a subsequence, the metrics converge smoothly to some limiting 
metric g = hi ^Ge, £ > 0 at least in a neighbourhood U of dCo , compare appendix A.2. 

Away from the points P± we can apply Lemma 3.7 to conclude that, after passing to a 
further subsequence, the maps converge on U* := U \ (Jj ± P± in the sense that 

Ui Woo weakly in H? oc (U*) and strongly in W l( ^{U*) for every p < 00 . 
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Furthermore, the uniform bounds on the energy imply that the maps iq converge to Uoo 
weakly in H 1 on all of U while the uniform Lr bounds on ||/i||i 2 give weak L 2 convergence 
to a limit foo again on all of U. Remark that here there is no need to specify with respect 
to which metrics gi the norms are computed as the metrics are uniformly equivalent. 

Furthermore the traces Ui\gc 0 converge uniformly to Moo|dC 0 thanks to the equicontinuity 
obtained from the Courant-Lebesgue Lemma, so that Uoo can be extended to an element 
of H^(C 0 ). 


We finally remark that the convergence of ( (v,oo,g, foo) implies that the 
differential inequality (3.12) is again satisfied for (uoo, g, foo) at least for variations sup¬ 
ported in U*, see also appendix A.l. 

The basic idea of the proof, working without modification only if the image of iq| gconB^ (x 0 ) 
happens to be the subset of a straight line, is now the following. 

Since Uoo is an element of H 1 , we can choose tq > 0 so that E(uoo, B 2ro (cco)) is far smaller 
than £2 and thus in particular far smaller than the energy of the maps iq on this ball. 
We would thus like to consider variations u e of which have the form tq +£ ■ X(uoo — iq), 

A a cut-off function supported on B 2 % ro {x o). Then, if we could insert v = ^| e= o u E as 
test-function into (3.3) we would get a contradiction since the first term would give a 
a negative contribution of roughly —£2 which could not be compensated by the second 
possibly positive term. 


Of course, having Plateau- rather than Dirichlct-boundary conditions, the maps rq + £ ■ 
A(woo — ut) are in general not in H p. To obtain an admissible variation we shall thus use 
the following lemma which is proved later on. 

Lemma 3.11. Let T be a regular closed Jordan curve in R" of class at least C 3 . Then 
there exist constants f = f(T) and C = C(r) such that for any point p £ T and any 
r € (0,f(T)) there exists a C 2 -diffeomorphism $ : 1" — > R" with <f ' = id outside of 
B 2 r{p) and 3>(p) =p which satisfies 

r 2 ||$ - id|| L oo + r||d$ - zd|U- + ||d 2 $||ioo < C (3.21) 

and which straightens out the curve T in a neighbourhood of p in the sense that 

$ : T 0 B r {p) -¥ {p} + T p T. 


Returning to the proof of Lemma 3.10 we let f € (0,f(r + )) be a fixed number that we 
determine later on. Then, as the traces Ui\gc 0 are equicontinuous we can choose ro > 0 
small enough so that 

Ui(B 9 2 ; o (x 0 ) n dC 0 ) C Br( Po ) (3.22) 

where po : = Ui(xo) = a + (0 o ) € R" is prescribed by the three-point-condition. 

Let now $ be the map given by Lemma 3.11 which straightens out the boundary curve on 
the ball Bf(po) C R” and let A € Cq° {B 2r {xo)) be a cut-off function which is identically 
1 on B^(x 0 ) C Cq and whose derivatives satisfy |Vg.Ai| < Cr^ fc , k = 0,1,2. Here the 
constant C is independent of i since the metrics gt converge smoothly near the boundary. 

We then define 

'1-4 := $ _1 o [$ o m + £ ■ Xi ■ (<h o Uoo — $ o ut)] 
and claim that, for £ > 0 sufficiently small, this is an admissible variation of iq, i.e. that 

u i e B^ ^(C 0 ). 
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As uf is clearly of class H 1 , and as u\ = Ui away from supp(Ai) it is enough to show that 
uf(xo) = p o and that uf \ac 0 r\B 3i (x 0 ) a wea kly monotone parametrisation of a subarc of 
r (namely of Ui^dCoDB^ (xq) )). The first claim is trivial since the three-point-condition 
is satisfied for all uf s and thus, since the traces converge uniformly, also for Uoo- The 
choice of ro and the properties of $ imply furthermore that the restriction of $ o m to 
8Cq (~1 Bf’ (a;o) gives a weakly monotone parametrisation of a segment in the tangent 
t Po ,r = {Po} + Tp 0 T- The same holds true also for $ o u^ and thus for any interpolation 
of these two maps so uf is indeed an element of *. Thus Wi := ^| £=0 uf € T+iJp * is 
an admissible test function for the differential inequality 

/ ( dui,dwi)dv gi + fi-Widv gi > 0. (3.23) 

J Co J Co 

This will lead to the desired contradiction for i sufficiently large provided the above 
construction is carried out on sufficiently small balls B®* and B f (po)- 

To begin with, we remark that 

Wi = ^| E=0 u\ = Aj • ((/<!'” 1 )(‘I>(?/,)) ■ ($(uoo) - $(u»)) (3.24) 

is supported in the small ball Bf" (a-'o) and bounded Hit?* ||x,°° < C2 by a constant de¬ 
pending only on T + and the bound M imposed on the L°° norms of the iq. 

We can thus bound the second term in (3.23) by 

I [ fi ■ Widv gi \ < C 2 ■ WMlubS. (*„)) ^ Cr o\\M\L2(c 0 ) < CAr 0 < f (3.25) 

J Co 

provided ro is chosen sufficiently small. 


The first term of (3.23) on the other hand can be bounded from above by 


J {dui, dwi) dv < J Xi ■ (d$ 1 )($(u i )) • d$(u 00) • du^ - dui 


, dui ) dv 


-I- C sup (|(d 2 $)(ui(a;))| • |<f>(u 00 (a;)) - $(uj(:c))|) • [ K\dui\ 2 dv 




'2r 0 


+ 


c 


L 


^0 Jsupp(dAi) 


|^(uoo) - $( u *)l • \dui\dv 


< ~ [l - ^ - C ■ w(r 0 )] ■ j \i-\dui \ 2 dv 

+ C • E(u oo, Bff. 0 ) + C • || Ui — Uoo 11 + CE{v,i , Bfro \ ^r' 0 ) 

(3.26) 

where all balls are to be taken with centre Xo , all integrals and norms are computed with 
respect to gi and where we set 

Wi(r 0 ) := sup (|(d 2 $)(ui(a;))| • |uoo(ar) - Uj(a:)|). 


x(zB 


2r 0 


Recall that m —*■ u^ in (U*) for every p < 00 and that the metrics converge. Thus 
the penultimate term in (3.26) tends to zero as i —> 00, and is in particular < j £ 2 for i 
large. Furthermore, for i large, the last term in (3.26) is bounded by C£'(u 00 , Bf (xo)) + 
j£ 2 < ^£2, where the last inequality holds provided ro is chosen sufficiently small. 

Given that \ J A, • |dui| 2 dv 9i > E{ui,B ®’) > e 2l we can thus estimate (for i large ) 

f (dui, dwi) dv 9i + [ fi ■ Widv 9i <-{}--CbJi(r 0 )) ■ [ Xi ■ \dui\ 2 dv 9i , (3.27) 
Jc 0 JCo 4 J 
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which leads to the desired contradiction to (3.23) provided we show that r 0 > 0 can be 
chosen so that 


CoJi(r 0 ) < — for i large. 


(3.28) 


To prove this last claim we recall that d 2 $ vanishes identically outside the ball B 2 f(p). 
This means that cc,; is obtained as supremum over a set on which the oscillation of the 
function Ui is a priori no more than 4f, for a number f that we can still reduce if needed. 
This aspect of the construction is crucial as we have no control on the behaviour of Ui 
near xq, so could in particular not hope for the oscillation of Wj over the full ball B^ to 
be uniformly small. 


For r > 0 sufficiently small and i large, we can in particular estimate 

Ccji(r 0 )<C- sup \ui-po\+C- sup | ^oo PO | 

B 2 r 0 ( ;,; o)nsupp(d 2 <I>Ou i ) B 2r 0 ( x o) 

< Cf + C • OSC Woo < — + C ■ OSC Woo- 

B lr a ( X o) 8 B 4r 0 ( x o) 


(3.29) 


We finally recall that w 00 satisfies (3.3) for the function /oo £ L 2 , the limiting metric g 
and for variations supported in U*. As such Corollary 3.8 and Remark 3.9 imply that 
Woo is continuous at least in a neighbourhood of dCo and thus in particular in the points 

Carrying out the above argument for a small enough radius ro, which might depend on 
Woo but is independent of i, we thus find that (3.28) indeed holds. □ 

It remains to prove 


Proof of Lemma 3.11. Let T be a C 3 closed Jordan curve, let po £ T C R n , let i POi r = 
Po + T po T be the tangent to T at po and let n : R” —> t POt r be the nearest point projection 
onto t POi r- 

Observe that for f > 0 chosen sufficiently small, in particular so that T n i?3f(p) is 
connected, this projection induces a C 2 bijection from T D B^f-ip) to a segment in t POt r- 
Furthermore, after possibly reducing r, we have that t P0 p fl B 2r (p) C 7r(r D B 3r (p)) for 
all r £ (0, f). 

We now consider if : t POy r fl B 2 f(p) —> R n defined by ip = id\t p T — (7r|rnB 3? (p)) 1 an d 
claim that given any number r > 0 

^ ^ r := id + X r ■ if o tt = id + X r ■ [V — (7r|rnB 3f .(p)) 1 ° 7r] (3.30) 

gives the desired diffeomorphism. Here X r £ CQ 3 (B 2r (p), R) is given by a cut-off function 
which is identically 1 on B r (p) and which satisfies the usual estimates of \D k X\ < Cr ~ k , 
k = 0,1, 2. Since Tr(TnB 2r (p)) C B 2r {p)C\t P} T C 7r(rn5 3 f (p)) the map is well defined 
for any radius r £ (0, f). 

To prove that T has the properties we asked for in Lemma 3.11 we first remark that 
7 x{p) = p and thus if{p) = 0; be. ^(p) = p. More generally, given any point x £ YC\B r (jp) 
we obtain that 

’F(x) = x + A r ■ [7r(ir) — x] = tt(x) 
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so U/ straightens the curve T D B r {p) to a line as described in the lemma. Since A r and 
thus also 'L — id is supported in B 2 r {p), it remains to show that the estimate (3.21) 
claimed in the lemma holds true with a constant independent of r. 

Since T is of class C 3 and since we project onto the tangent to T, an estimate of the form 
|7r(a:) — x\ < Cr 2 is valid for all x £ T D B 3r (p) (recall that T n B 3r (p) is connected). 

We then use that we can write any y £ i? 2 r(p) H t P0: r as y = ir(x) for an x £ T D B 3r (jp) 
to conclude that \if{y)\ = |7r(x) — x\ < Cr 2 . In particular 

ll'E' — id\\L°° < C sup \i/]ott\<C sup \if\ < Cr 2 

^2r{p) ^2r(p) CltpQ ,r 

holds true with a constant C depending only on T. We furthermore remark that since 
the derivative of the function if (which is defined only on a line) vanishes in the point p 
there exists a constant C (again depending only on T) so that ||cfo/>[|x,oo( t r ns 2 r(p)) — C r ■ 
We can thus bound 

Hd’I' - id\\ L °° < ||dA||ioo • \\ip o it\\ l °° + ||#||l~ • ||d7r|| L ~ 

< Cr- 1 -r 2 +C\\dip\\ L oc < Cr 

as well as 

< IM 2 a|| L co • ||^ o ttIU- + \\d\\\ L - ■ ||dv>||L- + c\\d 2 {ii> o ^)||i=o < c. 


□ 


An important consequence we can derive from our key Lemma 3.10 is 

Corollary 3.12. Let A, M < oo and let K be a compact subset of A4. Then to every 
e > 0 there exists a constant S > 0 such that the following holds true: 


Let u i and U 2 be such that (3.12) is satisfied for functions fi with ||/i||l 2 < A and metrics 
gi £ K. Suppose furthermore that E{ui,gf) < Eq and that ||uj||i=o < M, i = 1,2. Then 

if 


||wi - U 2 \\L^(C 0 ,g) < 8 


for some g £ K then also 


ll^i ^ u 2 \\m(c 0 , g ) < e. 


Proof. To begin with, we remark that all metrics in K are uniformly equivalent since K 
is compact. Thus it is sufficient to show the claim for norms || • ||^2 and || •II Hi that are 
computed with respect to some fixed g £ K. 

We then argue by contradiction. So assume there is a number e\ > 0 and triples 
(ui,<7i, fi)keN as well ( u 2.t92i /f )fceN so that all the assumptions of the lemma are sat¬ 
isfied (for each k) but for which 

ll«i - u^\\ L 2 (Co) 0 while || u\ - A\\ h hco) ^ (3-31) 

Then, after passing to a subsequence and using Proposition 3.5 and Lemmas 3.3 and 3.7, 
we find that locally on C* the maps u\ 2 converge to limits Mi .2 strongly in H 1 , where 
by construction these two limits must agree. 
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In particular, for any fixed number r > 0 we have 




(3.32) 


as k —> oo. 


We can then choose r > 0 so small that Lemma 3.10, combined with the equivalence of 
the metrics in K , implies that 

|K|K(b»(x 0 )) < ClKlL ? < K for all fc G N, i = 1,2, and x 0 G C 0 . 

Applied for the points and combined with (3.32) this contradicts (3.31). 


□ 


3.4 Convergence of the time-discretisation scheme 

Given any initial data (iio,So) G itp *(Co) x At we consider the approximate solutions 
of Teichmiillcr harmonic map flow ( Uj,gj) := (u hi ,g hj ), hj = 2~i , obtained by the time 
discretisation scheme described in section 3.1 We can analyse the maps Uj using the 
results of the previous section since Uj(t) can be seen as a stationary solution of 

J {duj(t),dw)dVg .( t ) + J Df 3 Uj(t) ■ wdVg^ t ) > 0 for w G T+ (t) Pp„ and t G [0,T], 

(3.33) 

where gj is piecewise constant so that gj(t) = gjit^ 1 ), t G [t*/Afc+i)- 

Based on the results of the previous sections we can pass to a subsequence, still denoted 
by (uj,gj), of approximate solutions which converge to a limiting curve of maps and 
metrics (it, g) as described below. 

To begin with, we claim that Uj converges uniformly in time with respect to L 2 in space, 
i.e. that 

sup || Uj(t) — u(t )\\ L 2 —> 0, for j — > oo (3.34) 

te[o,T] 

and that u G C°([0, T], P 2 (Co)). Here and in the following there is no need to specify 
with respect to which metric on Co the above convergence is to be understood as all the 
considered metrics are uniformly equivalent on the interval [0, T] we consider, compare 
(3.10). 

To prove this claim we let 1 K > Uj ( t ) be piecewise linear with Uj ( t ^) = Uj ( t ^) for every 
k. Then (3.9) gives uniform C t ' 2 L 2 estimates for Uj(t ), namely 

IKK) - KK)IK < £ WD^UjW^dt < (t 2 -H) 1/2 (^ \\Dt S Uj\\ 2 L3 dt) 1/2 

< (2Po) 1/2 K-H) 1/2 . 

Thus, after passing to a subsequence, Uj converges in C®L 2 to a limit it. Furthermore, 
again by (3.9), 

sup IKK — «jK|| < sup / \\Dt 3 Uj\\ L 2 dt < h 1 / 2 (2E 0 ) 1 ^ 2 —> 0 for j — > oo 

t k J t h k j 
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leading to (3.34). 


Remark that the uniform bounds on the energy furthermore imply that the limiting map 
is in L°°([0,T\, H 1 (M,go)) and that the spatial derivatives converge 

d,Uj —*■ du weakly in L 2 ([0,T] x Co). 

We recall that the traces of uj on dCo are uniformly equicontinuous, so we furthermore 
have that 

Uj(t)\gc 0 —> u(t)\dc 0 uniformly on dCo for all t £ [0,T], 
where we stress that we do not claim that the rate of this convergence is uniform in time. 

Additionally, the energy inequality (3.9) gives uniform L 2 (Cq x [0,T — hj]) estimates 
for the difference quotients D^ 3 Uj. Consequently, u is weakly differentiable in time on 
Co x [0,T] with D^ 3 Uj —*■ d t u in L 2 (Cq x [0,T]). 

For the metric component gj we can apply Lemma 3.2 to get uniform C 0,1 estimates in 
time with respect to any C k metric in space since the L 1 norm of the Hopf-differential 
is bounded in terms of the (non-increasing) energy. 

We can thus get convergence of gj —*■ g in C o, “([0, T], C fc (Co)), a < 1, with the limiting 
curve being again of class C t u ’ C k and thus in particular differentiable in time for almost 
every t £ [0, T\. 

We shall now prove that the limit (u, g) obtained in this way gives the desired solution 
of Teichmtiller harmonic map flow, namely that 

Proposition 3.13. Let (u hj , g hj ) be a sequence of approximate solutions to a fixed initial 
data (uo,go) converging as described above to some limiting curve ( u,g ) as hj 0. Then 
the limit (u, g) is a stationary weak solution of Teichmuller harmonic map flow which 
also satisfies the energy-inequality (for a.e. t\ < t^)- 

We remark that while g is clearly again an admissible curve, we need to prove that its 
derivative is actually given by the projection of the Hopf-differential of the limit. As the 
projection operator is non-local, for this part of the proof the key lemma 3.10 and its 
Corollary 3.12 are crucial to get strong H 1 convergence for the map u and thus strong 
L 1 convergence of the Hopf-differential on all of Co, in particular also near the points 
Pp where Proposition 3.5 does not apply. 

Conversely, the analysis of the map component can be carried out very similarly to 
the work of Duzaar and Scheven [6] and is indeed less involved then the corresponding 
arguments since the metric is well controlled and since our equation for the map is linear. 


Proof of Proposition 3.13. We first infer from the energy inequality (3.9) that for any 
A < oo the set of times 

A ^ := {t £ [0,T] so that \\Dt 3 Uj(t)\\ L 2 < A} 

has measure £ X (A^) >T— ^r, so in particular £ 1 (A a ) >T— also for 

OO OO 

A a = limsup A a = P| A a . 

n=l j=n 
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Recall that the maps Uj(t) satisfy (3.12) for / = D^ 3 Uj(t) so it is precisely bounds of the 
form ||nf J w J (t)|| i 2 < A that are required in order to be able to apply the results derived 
in the previous section. 

We begin by analysing the metric component. To prove that g indeed solves (2.7) we 
show that it agrees with the solution g(t) € of the initial value problem 

dtg = P§{Re((&(u,g)), g(0)=g o . 

We will prove this claim based on Corollary 3.12. So let K be the set of metrics hi ^Gt 
satisfying (3.10) as well as (3.11) and let M be a bound on the L°° norm of the initial 
map which therefore also serves as bound for compare Lemma 3.1. 

Then given any numbers e > 0 and A < oo we let S > 0 be the number given by Corollary 
3.12 and select jo so that 

sup || Uj(t) - u k (t )\\ L 2 < S for j,k > j 0 - 

te[o,T] 

This implies that 

\\u k (t) - Uj(t)\\ H i( Co ) < £ for all t G A a n A A and j, k > jo 
which in turn yields the same bound for \\u(t) — Uj(t) ||hi(c 0 ) f° r t S A A n A A . 

For times in A A fl A A the difference of the Hopf-differentials is thus controlled by 
\\Reg(<f>(u,g){t)) - Re gj (^(uj,gj)(t))\\ L i < C ■ ||(ft ~ g)(t)\\c° ' E 0 

+ C-E$.\\(u j -u)(t)\\HHC 0 ) ( 3 ' 36 ) 

< C- \\(9j ~ 9)(t)\\c° +C-e, 

c.f. (3.7). 

Here and in the following constants C may depend on Eq, K and T but not on e or j 
unless indicated otherwise. 

Based on Lemma 3.2 we can thus conclude that 

\\dt(g - g 3 ){t)\\ ck < C ■ \\(g - gj)(t)\\ C k +Ce 
for any t € A A fl A A and for j > jo (A, e). 

On the other hand, we can always bound the norms of dtgj and dtg by C • Eo and we shall 
use these trivial bounds on the set [0, T] \ ( A A . fl A A ) on which we cannot apply any of the 
results of the previous section since we lack the necessary control on the inhomogeneity 
of (3.12). 

Combining these two cases, we obtain that for almost every t 

\\dt(9- gj)(t)\\c>° < C- \\{g -gj)(t)\\ c * +h A (t) 

where 

h A = C ■ e + Ci?oX[o,T]\(AAnA A )- 

Based on Gronvall’s Lemma, we can thus conclude that for any t £ [0, T] and j > jo(e, A) 
ll(fl-flj)(*)llc* <e CT ■ \h A \dt<Cs+ 
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Choosing A —>• oo and e — > 0 and corresponding values of jo(£, A) —> oo yields the claim 
that gi —> g uniformly and thus that g = g is indeed the solution of (2.7). 

We now turn to the analysis of the map component where we follow largely the arguments 
of [6], 

To begin with, we observe that for almost every time t £ [0, T] there exists a number 
A < oo such that t £ A A . Choosing a subsequence along which \\Dt 3 Uj(t)\\ L 2/c 0 ) —»• 
liminfj-joo WD^ 3 Uj{t)\\L^tc 0 ) < A we conclude that 

Uj {t) -> u(t) 

converges not only strongly in L 2 , but thanks to Corollary 3.12 indeed strongly in H 1 on 
all of Co and, thanks to Lemma 3.7, also weakly in H 2 oc (C*) where C* := Co \ 

We stress that the choice of this subsequence is allowed to depend on the time t we are 
considering. 

We furthermore remark that combining the uniform P 2 -estimates for Uj valid on sub¬ 
sets O CC C* with the uniform convergence of the metrics gj —> g yields that also 
A gj (t) u j{t) — A g (t)u(t) weakly in L 2 (Vl) for each such fh Using the L 2 bound on 
A g .( t ) u j(t) = D^Uj valid for t £ A A we can thus conclude that 

l|A s ( t )w(f)|| L 2 (n) < lim \\A g {t )Uj(t)\\ L 2 ( n) < A for each Q CC C*. 


Passing to the limit in both the Euler-Lagrange-equation and the stationarity condition, 
cf. appendix A.l, thus yields that u(t) is a stationary solution of (3.12) for a function 
f(t) = A g ( t )u(t) whose L 2 (Co,g)-norm is again bounded by A, or indeed more precisely 
by liminf-^oo \\D^ Uj{t)\\ L ^ Co , gj y 

Repeating this argument for a.e. t £ [0, T] we thus obtain a function / : [0, T\ x Co which 
must have bounded L 2 (Co x [0,T])-norm since 

ll/lll 2 (C 0 x[o,T]) < [ lim inf \\D^ Uj(t)\\\ 2{Co) dt < liminf ||D t ,tj M i ||| 2(CoX[0)T]) < 2 E 0 
Jo J —* 00 J—^oo 

where the last inequality follows from (3.9). 


We now wish to show that / agrees with the time derivative of u. To this end, proceeding 
as in [6], we set 


uf(t) 


u(t) if f £ B a := [0, T] \ A A - a 
Uj (t) if t £ A A and fj (t) — 


fit) if t £ Bf 

Dt j u j (t) if t£A A 


in order to obtain a new sequence of pairs satisfying (3.12) for gj(t) but for which the 
estimate ||/j A (t)||L2(Co) < A is now satisfied for every j and every t £ [0,T]. 

We first claim that for any sequence A j —> oo 

ff ' 3 d t u weakly in L 2 (Cq x [0,T]) for j — > oo, 


or, as we know that D^Uj —- dtu, equivalently f^ : — D^Uj —^ 0. Indeed, given any 
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function cj) £ L 2 (Cq x [0,T]) we have that 


' [0,T] xCo 


</>■ ( D t iu i ~ fj J )dv g dt < ll^ll i2(B A i xCo) • (II//IIl2([0,T]xC o ) + \\ D t j Uj\\L^[0,T]xGo)) 


<cE 2 m T2f , 


L 2 (Br 3 xC 0 ) 


(3.37) 


which tends to zero as Aj —> oo since the measure C{B - 3 ) —> 0. 


At the same time we claim that du■ 3 converges not just weakly, which would be evident 
from an argument just as carried out above, but indeed strongly in L 2 ([ 0, T] x Co) to du. 
Indeed, let us first consider du — du A for a fixed number A. Since this difference vanishes 
on and since we can apply Corollary 3.12 to u(t) and uf{t) for j large enough (so 
that these maps are L 2 close) we obtain that 

II du - duj\\ L 2 iCoX[0tT ]) 0 for every fixed A. (3.38) 

But the set of times AA a c B^ 3 U B^, on which du A and du^ 3 do not agree, has 
measure no more than C(A~ 2 + A~ 2 ). Combined with the uniform bound on the energy 
this means that (3.38) suffices to conclude that indeed || du — du^ 3 ||l 2 (c 0 x[o,t]) —> 0. 


Thanks to the strong H 1 convergence we can furthermore approximate each test function 
w £ L 2 ([0,T],T+77p *(Co)) by elements Wi £ L 2 ([0, T], T+iJp *(Co)) in the sense that 
\\dw — dwi\\i J 2 ^ CoX ^ O T ^ —> 0, compare appendix A.l. We thus conclude that 

I du ■ dw + dtu ■ wdv a dt > 0 
Co 

for all such w. As we have already shown that g satisfies (2.7), we thus obtain that (u,g) 
is indeed a weak solution of Teichmiiller harmonic map flow. 



Knowing that du^ 3 converges strongly and not just weakly furthermore implies conver¬ 
gence of the Hopf-differentials $(uj,gj) to <F(u, <7) in L 1 , allowing us to pass to the limit 
in the stationarity condition to conclude that (2.11) holds true for almost every time t. 


It remains to show that the energy inequality holds true (for almost every pair of times 
t\, t2-) 

So let t\ £ [0,T] be any time which is contained in one of the sets A A , A < oo, i.e. for 
which there exists a sequence of hj so that D^ 3 Uj{t\) is bounded in L 2 . For this sub¬ 
sequence of Uj(t i), Corollary 3.12 implies strong H 1 convergence and thus in particular 
that E((uj,gj)(t!)) E((u,g)(t i)). 


We then recall that Uj satisfies the energy inequality 

~ E(uj,gj)(t) > / \\dtgj\\L 2 (c 0 ,gj)dt + 2 / ll^* u jll L 2 (M,gj(t))dt 

J t\ J t\ 

for all times t £ [0, T], t > t\ and Uj(t) converges at least weakly in H 1 to u(t) (a further 
strongly convergent subsequence could be found for almost every t but is not needed). 
We thus obtain that for every t>t\ 

E(u,g)(h)-E(u,g)(t) > ( \\dtg\\ 2 L 2 (Co ,g)dt + \ [ \\d t u\\ 2 L 2 {c ^ g) dt. 

Jtx j tl 

□ 
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4 


Long time existence 


4.1 A priori estimates for the metric component 

Before we can analyse admissible curves of metrics in more detail we finally need to 
decide how to select the family of diffeomorphisms hi ht j, which we use to compensate for 
the lost degrees of freedom of the three-point-condition. 

Rather than just writing down a possible family, we shall first describe which properties 
we require in the present context of flowing to minimal surfaces. We will then later give 
an example of such a family but do not claim that this choice is in any way unique. 

To begin with, in order to obtain solutions of the flow that exist for all times we need 
the following L 2 -completeness property: 

Lemma 4.1. Let ( hi,he the family of diffeomorphisms defined in (4.2). Assume that 
9{t) = ^b(t) t £ [O,! 1 ), is such that the diffeomorphisms hb ,0 become singular as 

t T, i.e. so that (at least) one of the values | 6 ± | —»■ 1 or |</) ± | —> oo as t —»• T. Then 

\\dtg\\L*(c 0 ,g)dt = oo. 

A further requirement we want to impose in preparation for the asymptotic analysis 
carried out later in section 5 is 

Lemma 4.2. Let (hb,f) be the family of diffeomorphisms defined in (4.2) and let y(&, <Jf) 
the space of generating vectorfields of hb^- Then 

T(TCo) = T(TC'o)* © /i^x(M) 

for all (b,<j>) £ f lh '■= D i(0) 2 x R 2 . 

Here and in the following y( 6 , <j>) C T(TCq) is the 6 dimensional vectorspace spanned by 
the vectorfields generating the diffeomorphisms hb , 0 , i.e. by Y ( p± (b, <p) characterised by 

d^hb^is, 9) = F 0 ± ( 6 , f)(h b ^{s, 9)), (4.1) 

together with the vectorfields Y Re ^ b ±)(b,(j)) and Y Im ( b ±){b, <j>) defined by the analogue of 
(4.1) or, if b ± 0 , equivalently together with the vectorfields Yj h ±|( 6 , <j>), Y Arg tb±)(b,<t>) 
corresponding to variations of the absolute value respectively the argument of b ± . 

The final property we shall ask of the diffeomorphisms h bt< j> is that their support is disjoint 
from the middle geodesic. This will have the advantage that the modification by these 
diffeomorphisms does not interfere with the analysis of a possible collapse of the central 
geodesic, cf. Lemma 4.4. It will furthermore prove to be useful to choose the hb,</> so that 
the support of the induced variations of the metrics with respect to the parameters 
on the one hand and b ± on the other hand are disjoint. 



4.1.1 Choice of diffeomorphisms 

A simple way of assuring that our diffeomorphisms satisfy Lemma 4.2 is to choose them 
as restrictions of Mobiustransforms on the boundary of Cq. 
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Given numbers 0 ± S M and 6 ± eC with | 6 ± | < 1 we consider the functions f b ± j( />± '■ R —> 
R which are induced by the Mobiustransforms M b ±^± , i.e. chosen so that /o,o = id and 
e *f b ±' 4 ,±( 0 ) — M b ±^±{e ie ), where 


M b Jz ) := Z+ - b , for b e £>i(0) C C, 0 € R, z e £>i(0) C C. 

’ 1 + 0-2 

We then extend the maps induced by f b ± :< j,± on OCq to a suitable diffeomorphim h b ,0 on 
the whole cylinder. Namely, we choose Ai .2 as smooth cut-off functions such that Ai = 0 
on [—1, |] with Ai = 1 on [|, 1] while A 2 = 0 on [—1, A] with A 2 = 1 on [|, 1], We then 
define h b ^ : Go -+ Co, 6 = (b~,b + ), (f) = (</>“, (j) + ) through 

hbA*, 0) = (a, Ai(s) • / 6+ ,*+ ( 0 ) + (1 - A^)) • (t9 + A 2 (s) • 0+)) (4.2) 

if s > 0 respectively by the analogue formula, replacing (& + , </> + ) with ( 6 _ , <(> _ ) and s by 
—s, if s < 0 . 


Since f b ,<j>{Q) = fb,o(d)+cj) this formula reduces to h b ^(s,d) = (s, 0+\i-(f b {0) —9)+\ 2 ■(/>+) 
where we write for short f b for f b ,(j>= o- 

In order to show that this family of diffeomorphisms satisfies Lemma 4.1 we observe that 
a change of one of the parameters, say of |fo + |, induces a change of the metric of 

l^\( h b,4> G ) = K,<p L Y ]b+] G = L K ^ b+] g 

for g = h b ^ G and that the resulting Lie derivatives of the collar metrics G satisfy the 
following estimates 

Lemma 4.3. Let ( h b>t f >) be the family of diffeomorphisms defined in (4.2), let Y| b ±|, 
^ / Arg(h±) an d be its generating vectorfields and let (Gf) be the family of metrics 
defined in Lemma 2.4 for some fixed number 77 > 0. Then to any number Lq < 00 
there exist constants C\ y 2,3,4 € R + (depending only on Lq and g) such that the following 
estimates hold true for any metric G = Gn with £ < Lo and any ( b , (f>) (where we assume 
that b + 0 for the first estimate) 

||-£'Y | 6 +| G|| L 2 ( CoiG) > ^ _ |k + 1 - C 2 and C 3 < \\L Y<t>+ G||z, 2 ( GoiG) < G 4 (4.3) 
Furthermore 

Ly 1 G, Ly , G and Ly , G 

J |6+| ’ 'Ar s (li+) 1 <t>+ 

are L 2 (Gq, G)- orthogonal to each other. 


The claims made above for variations with respect to </> + and b + are of course valid also 
for variations with respect to b~ and <f>~ and from the construction it is evident that 
variations with respect to (</> + , b + ) on the one hand and {(j>~ , b~) on the other hand have 
disjoint support so result in Lie-derivatives that are trivially orthogonal. 

With regards to the proof of this lemma, we observe that the orthogonality of Ly. + G to 
the variations with respect to b + follows since Y^+ is given by the Killing field on the 
support of Y| b +| and Y Arg ( b +y 

The orthogonality of Ly^ + ^G and Ly A {b+) G on the other hand will follows from the 
different symmetry properties of these two tensors 
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The proof of this last part and of the estimates claimed in the lemma is not difficult 
though a bit technical so we include it in the appendix A.3. 

As a consequence of Lemma 4.3 we can now prove Lemma 4.1 for this particular choice 
of diffeomorphism 

Proof of Lemma 4-1- Let g(-) = ^b(.) <j,(.)Gt{-) be an admissible curve of metrics with 

L Lig(-) = fo \\ d tg\\L*(c 0 ,g) dt < °°- 

We first recall that Re(3i{g)) is orthogonal to {Lxg} so that both ||9 t G|| i 2(c 0 j G) < 
\\dtg(-)\\mc 0 ,g) and ||^| £ =o/i6(. +£ )^(. +£ )G(-)||i2 (Coi9) < \\dtg(-)\\ L ^Co,g) must have finite 
integral over [0,T). 

On the one hand, this implies that iff) is bounded from above by a constant L depending 
only on the initial metric and L^gf), compare (A.9). 

Using the orthogonality of Ly + G to the variations generated by a change of any of the 
other parameters, as well as estimate (4.3), we know furthermore that 

II^IU’(Co, 9) > \&<t> + \-\\4? h U G \\ l* (CoKa g) = \i^\-\\L Y , + G\\ LW) > C 3 -\icf+\ 

where C 3 depends only on the upper bound on £ obtained above. This implies that (f > + , 
and by the same argument also (j>~, remains bounded. 

So consider instead the behaviour of b ±, say of b + . The orthogonality relations of Lemma 
4.3 combined with (4.3) imply 


> iiAi^I ■ Lr^p\\W„,G) > [tAv ^ a] 2 iii4+n 2 . (4.4) 

In particular, for |6 + | sufficiently close to 1, an estimate of the form 

\£t l°g(l “ l 6+ DI ^ c \\ d tg\\L^M, g ) 

holds true which prevents b + from reaching dD i(0) if the curve g has finite L 2 length. □ 

We remark that the Teichmiiller space of the cylinder equipped with the metric that 
results from representing conformal structures by hyperbolic metrics f*Gj as described 
in Lemmas 2.1, 2.2 and 2.3 is not complete. Indeed, as explained in appendix A.2, for 
general curves in A4 and £ small we can only bound 

lfl<C'-||5 t5 |U 2 -£V2 

so that the possibility that l —> 0 is not excluded for curves of finite length. 

Nonetheless, for Teichmiiller harmonic map flow a degeneration of the metric in finite 
time is excluded since we can prove 

Lemma 4.4. To any numbers t\ > 0 and M, T, Eq < oo there exist constant C < oo 
and Eq > 0 such that the following holds true. Let (itojffo) 6 -^r*(Co) x A 4 be any 
initial data so that E(uo,go) < Eq, UmoIU 00 < M and inj(Go,<?o) > 2^i and let ( u,g ) 
be the corresponding stationary weak solution of Teichmuller harmonic map flow whose 
existence on some interval [0,71) is assured by Proposition 3.13. 
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(4.5) 


Then the weighted energy is bounded by 

I(t) := [ e(u{t), g{t))p~ 2 [t)dv g{ t) < C, 

J Co 


and the injectivity radius by 

inj(Co, <7) > e 0 

for every t £ [0, min(T, Tj)). Here e(u, g) = ||ekt| s is the energy density while p(t)(x, 9) = 
Pt(t)i s e(t){ x )) Is the conformal factor of the hyperbolic collar. 


This result is essentially a consequence of results proven in [15] for Teichmiiller harmonic 
map flow from closed surfaces into non-positively curved targets because the action of the 
diffeomorphisms hb^ does not affect the region near the central geodesic. We also recall 
that while the metric component g is in general not smooth, it is Lipschitz continuous in 
time with respect to any metric in space. So while u might not be smooth in the interior 
of Cq it also satisfies such Cf’ 1 C^ bounds, at least away from time t = 0 which is enough 
to apply the arguments of [15] on almost every time-slice. 


Proof. We first explain why a bound on the weighted energy I results in a bound on the 
injectivity radius. We recall that the evolution of g(t) = h b,*m G ®> G W = G ^)> 
splits L 2 -orthogonally into the projection of the Hopf-differential onto the subspace 
{Lh^xg, X £ x(&!0)} and into the projection onto Re(H{g)) and thus that d t G = 
Re(P HI ' G \<&)). For the cylinder the space TL{G) consists only of tensors that can be 
written as 

do • dz 2 , ao £ K 

with respect to collar coordinates z = s + i9, (s, 6) £ [— Y{£), Y{£)] x S' 1 , so 

(46) 

We furthermore recall that if dtg = aodz 2 , then the length of the central geodesic evolves 
according to ^ compare (A. 8 ). 

For small values of i, say £ £ (0,£o)i the norm ||Gk 2 ||| 2 is given by (A. 6 ) which, once 
combined with (4.6) and (A. 8 ) and (A.4), implies that 

l4 log ^ + 7^-3 f (K| 2 - \ug\ 2 )p~ 2 dsd9\ < C£\\<f>(u,g)\\ L i, 
dt 16tH J c{t) 

C(£) = [— Y{£),Y{£)\ x S 1 , compare also section 5 of [15]. 

In particular log(£)| < C£+£I(t) is bounded if I(t) is bounded, resulting in the desired 
lower bound on £ = 2inj(Co, Gf). 

For the proof of (4.5) we can use results derived in sections 3 and 5 of [15]. Namely, the 
results of [15, section 3], in particular Proposition 3.6, give angular energy estimates for 
maps from hyperbolic collars into compact non-positively curved targets. Since we know 
that ||u||£«. < M these results apply without change also to the present situation. 

As in section 5 of [15] we consider a cut-off version of the weighted energy given by 

l(t) := [ e{u{t),g(t))p- 2 (t)ip(p(t))dv gW , (4.7) 

J cm) 
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where ip £ Co°([0, 2(5), [0,1]) is a cut-off function with p = 1 on [0, <5], and where 5 > 0 
can be chosen to be any fixed number. 


We remark that we can choose S sufficiently small, so that the diffeomorphism hb,$ agree 
with the identity on the support of <p o p, compare (A.3) and the subsequent comments. 

For such a choice of S we conclude that the evolution of the metric reduces to dtg = 
Re(c(t)dz 2 ) on the relevant region, i.e. on the support of p o p. 

Consequently the Bochner formula for the energy density given in Lemma 5.2 of [15] and 
the evolution equation for the conformal factor described in Lemma 5.4 of [15] apply 
without change and could indeed be further simplified as dtg evolves not just by any 
holomorphic quadratic differential but by cgdz 2 . 


Then arguing precisely as in the proof of Lemma 5.1 in [15] we obtain that 


dt 


log(l+I) 


< (7 (^1 + ||AgM||^ 2 ( Co S )^ < C (l + ||dtw||^ 2 ( CoiS )) (4.8) 


with C depending only on M, the initial energy and the choice of 5. Thus I and 
consequently also I < X + CgEo is bounded uniformly on every compact time interval as 
claimed in the lemma. □ 


From Lemma 4.4 we thus conclude that for arbitrary initial data (uo,go) £ x A4_i 
solutions to Teichmiiller harmonic map flow from the cylinder indeed extist for all times 
as claimed in Theorem 2.6. 


5 Asymptotics of global solutions 


We now turn to the proof of the second main result of the paper, the asymptotic conver¬ 
gence for the global weak solutions whose existence we have just proven. In the present 
work we analyse the asymptotics in case that the three-point-condition does not degen¬ 
erate as t —> oo, i.e. for solutions for which the parameters b ± remain bounded away 
from dDi (at least for a subsequence tj —> oo). The remaining case of the asymptotics 
will be analysed in future work. 

So let (u, g ) be a global stationary weak solution of Teichmiiller harmonic map flow 
which satisfies the energy inequality. We then choose ti —> oo such that the stationarity 
condition is satisfied for the times t,; and so that 

II A(ti) u (^)ll^ 2 (Co,s(ti)) °> ( 5 - 1 ) 

\\pV(Re^{u,g){U))\\ L . { c 0 Mu)) 0 (5.2) 

and |5 ± (ti)| -t* 1 as i — > oo. 

We can thus pass to a subsequence to achieve that 

bf = b ± (t i ) G £>i(0) C C and <j>f := 4>f - nf ■ 2w (5.3) 

as i —> oo where rii = j. 

We then pull-back the map and metric by the diffeomorphisms ft := ■ Remark 

that since /; = h^. ^ ° h b _ the resulting metrics gt := f*g(ti) are simply given by 
9i = 
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We furthermore recall that /,; agrees with the identity in a neighbourhood of the boundary 
so that the pulled-back maps Ui = u(U) o /) still satisfy the three-point-condition, i.e. 
are again elements of H p t (Co). 

Remark that (5.1) and (5.2) are satisfied also for (iii,gi) and both the differential in¬ 
equality 

J (dui,dw) dv 9i + J A gi m ■ w dv 9i > 0 for all w £ T+Pp* (5.4) 

and the stationarity equation 


Re($(ui, gi)) ■ Lxgt + A gi m ■ dui(X) dv g = 0 for all X G r*(TCo) (5.5) 


hold true. 

To prove convergence of ( Ui,gi ) to a critical point of area as described in Theorem 2.7 
we now distinguish between the non-degenerate case, i(ti) 0, in which wc will obtain 
a (branched) minimal immersion parametrised over a cylinder, and the degenerate case 
£(ti) —> 0 in which the surface splits into two minimal discs. 

We begin with 


Proof of Theorem 2.7 part (i): The non-degenerate case. After possibly passing to a fur¬ 
ther subsequence we can assume that ti = i{tf) —> £oo > 0 which implies that the metrics 
converge g t ->• g^ = smoothly on C 0 - 

Furthermore, as Ui is a solution of (5.4) for which ||A 9i Ui || i 2 is bounded, we can apply 
the if 2 -estimates of Lemma 3.7 away from P- as well as the H l estimates of Lemma 
3.10 and Corollary 3.12 on the whole of Cq. We conclude that a subsequence of the iq 
converges to a limit Uoo £ Hf oc {C*) D H 1 (Cq) where the obtained convergence is weak 
Hf oc and strong convergence on C* := Co \ (J Pf as well as strong H 1 convergence 
on all of Co- Furthermore, Corollary 3.8 implies that the maps iq are equicontinuous 
near the boundary, and thus by the H 2 estimates on all of Co, so that the ut converge 
uniformly on Co- In particular, u Q0 £ C°(Co). 

The above convergence implies not only that 

A 9oo Uoo = 0 on C 0 (5.6) 

and consequently that the Hopf-differential of the limit is holomorphic, but furthermore 
that the Hopf-differentials $(«i,<7i) —► $(«oo,3oo) converge in L 1 on the whole cylinder 
C 0 . 

From Lemma 3.2 we thus obtain that 

P sL( $ ( u °°>5oo)) = lim P 12 ($(u i ,3 i )) =0. (5.7) 

^ 2—>-00 a 


On the one hand, this implies that 


J Pe(<F( 

^OO 1 9 OO )) • L vg o 


,dv goo = 0 


(5.8) 


holds true for the vectorfields Y £ h* - x(boo,<t>oo ) generating the diffeomorphisms 

"OO 5000 
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On the other hand, the convergence of the Hopf-differential allows us to pass to the limit 
in the stationarity condition to conclude that (5.8) holds true also for all vectorfields 
X £ r(TC'o).. Thus, by Lemma 4.2, we find that (5.8) is indeed true for any smooth 
vectorfield on Co which is tangential to dCo on 8Cq. 

We now show that this forces 4*00 = $(«<», ffoo) to be of the form cdz 2 for some c £ R. 

Remark that if $00 were smooth (or even just W 1 ' 1 ) upto the boundary, we could directly 
combine (5.8) with Stokes theorem to conclude that $00 is real on the boundary and then 
to conclude that $00 = cdz 2 , c £ R. 

However, while $00 is holomorphic and thus smooth in the interior as well as in W 1,p , 
p < 2 in a neighbourhood of general boundary points, near the points P- we know a 
priori only that $00 is in L 1 . Thus $00 could have a pole at such a point and we need 
to proceed with more care. 

Given any fixed X £ T(TCq) we use that (5.8) implies that 


Lxg 00 Pe(-F 

OO ) d v g 


0 as € —0 


J [—l+e,l—e] xS 1 

and we initially work on such subcylinders where $00 is smooth. 


(5.9) 


Recall that Lxg can be identified with —5*X, where S* is the L 2 -adjoint of the divergence 
operator and that the real part of a holomorphic quadratic differential is divergence free. 

So, switching to collar coordinates (s,8) £ [—Toe, Wo] x S 1 , Yoo = Y(i 00 ), and applying 
Stokes theorem to (5.9) yields 

| f Re($ 00 )(-^,X)p~ 2 dd - f Re(^ 00 )(-§ i ,X)p~ 2 d6\ —s> 0 as i ->■ 0 

Jn' -eixS 1 it-y^+etxS 1 


'{T^-efxS 1 
where p = pt ^ ( s ). 


(5.10) 


Away from the boundary of [-Y.Y] x S 1 we now represent $ by its Fourier expansion 
•Foo = J2nez( Un + ib n )e ns e nie , a n ,b n £ K and apply (5.10) for vectorfields of the form 
X = A ± (s)cos (md)--jpg andX = A ± (s)sin(m 6 1 )-^, m £ N, where A ± are cut-off functions 
that are identically one in a neighbourhood of ±1 and that vanish say on {±s < ^}. 

Passing to the limit e —>• 0 in (5.10) yields 


h p 
u m c - 


mYoc 


= b„ m e~ mY °° and b m e~ mY ~ = b m e mY ~ 


as well as 




= —a_ ? 


,e _rn Y “ and a, 




so that all Fourier coefficients except for Co = ao + need to be zero. Of course, testing 
with X = X ± ■ furthermore gives that 60 = 0 and thus that (Foo = a^dz 2 is indeed an 
element of R(Cq). 

But (5.7) also implies that the projection of $00 onto R{goo) = {cdz 2 ,c £ R} vanishes 
so $00 must vanish meaning that u 0c must be (weakly) conformal. Thus Uoo is a weakly 
conformal and harmonic map which spans T and can thus in particular not be constant 
so must be a (possibly branched) minimal immersion [7]. □ 


Proof of Theorem 2.7 part (ii): The degenerate case: Let (m, gi) be as above and assume 
now that ii —> 0. We let C + = (0,1] x S 1 and C~ = [—1,0) x S 1 and observe that the 
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subcylinders (C' ± ,gj) are isometric to 

(M) x S l ,(${Yi - s) ■ (ds 2 + d9 2 )), Yi = Y(£i) 

with an isometry given by ff; : (x, 6) i-»- (Yj =p (a;), 0). We remark that pt(Y (t) — s) —> 
2ir l+ v as l —> 0 locally smoothly on [0, oo) x S' 1 . At the same time converges locally 

to the diffeomorphism : C ± —> [0, oo) x S 1 given by f^(x, 9) = (|^ • tan (f =F — 
2nr], 9). 

Thus the metrics gt, converge smoothly locally to a metric goo that is isometric to the 
hyperbolic cusp 

([0, oo) x S 1 , Pq(s) • (ds 2 + d9 2 )) 

described in the theorem. 

At the same time, we get subconvergence for the maps iq = u(U) o ho t 2 nni —l Uoo as 
described in the theorem since the bounds on \bf\ allow us to apply the i7 2 -estimates of 
Lemma 3.3 and Lemma 3.7 as well as the H 1 estimate of Lemma 3.10 and the equicon- 
tinuity result of Corollary 3.8 on every compact subset of C ± , see also Remarks 3.6 and 
3.9. 


The above convergence of the maps and metrics implies in particular that Auoo = 0 
and thus that $00 = ^(uoo, goo) is holomorphic on C±. We then observe that the local 
convergence of the map and metric on C± allows us to pass to the limit in the stationarity 
condition to conclude that 



Txffoo • i?e(i>oo)du Soo 


= 0 


(5.11) 


provided we only consider vectorfields X £ r(TCo)* whose support is contained in one 
of the subcylinders C ± . 


We recall that also the support of the vectorfields Y^±. Y Re ^ b ±) and Xrm( 6 ±) generating 
the diffeomorphisms h b ^ is contained in and that the corresponding projection of $ 
tends to zero, compare (5.2). So local strong convergence of Ui in H 1 and consequently 
of in L 1 implies that (5.11) holds true also for these particular vectorfields, and thus, 
by Lemma 4.2, indeed for arbitrary vectorfields X £ T(TCo) whose support is contained 
in either C + or C~. 


As the Fourier expansion of &(uoo, goo) on (C , g^) ~ [0,oo) x S ' 1 cannot have any 
exponentially growing terms (since ||4 > 00 ||ii < CEq < 00 ) we can then argue as in the 
previous proof to conclude that in collar coordinates (s,9) £ [0,oo) x S 1 

‘hoo = c ± (ds + id9) 2 for some c ± £ R. 


Pulling $00 back to the punctured disc D* = D i(0) \ {0} through a conformal diffeomor¬ 
phism / : ( D*,g euc i ) —> ([0,oo) x S 1 ,ds 2 + d9 2 ) we thus find that the Hopf-differential 
of Uoo is represented by c^z^dz 2 , £ R for z £ D*. 

But the limiting map Uoo can be seen as a harmonic map from the punctured disc 
(D *, g euc i) whose energy is finite (since the energy is conformally invariant) which implies 
that Uoo can be continued smoothly across the puncture, compare [16]. Thus <I > 00 must 
be smooth on all of -Di(O) and must thus vanish identically. 

This proves that the maps = Uoo|c± extend to weakly conformal harmonic maps 
from the disc and thus give two (possibly branched) minimal immersions with each of 
them spanning one of the boundary curves T . 
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□ 


A Appendix 

A.l Courant-Lebesgue Lemma and properties of Hy. *(C 0 ) 

Throughout the paper we made use of the Courant-Lebesgue Lemma of which we use 
the following version, see e.g. [9, Lemma 3.1.1] or [17, Lemma 4.4] 

Lemma A.l. Let D r ( 0) + = {x £ R 2 : |x| < r,x i > 0} and let u £ H 1 (D r ( 0) + ,R ra ) 
be any map that has energy E(u,g euc i) < Eq, Eq any fixed number. Then for any 
<5 £ (0,min(r, ^)) there exists p £ (S,\/6) so that u\qd ( 0 )+ is absolutely continuous and 
so that the estimate 

\u{x) - u(y)\ < C- | log((5)r 1/2 , for all x,y £ dDj := {y : \y\ = p, yi > 0} 
holds true with a constant C that depends only on Eq. 

We use in particular the following consequence for maps satisfying the three-point- 
condition 

Corollary A. 2. Let Ui £ Hf * (Co) be a sequence of maps that have uniformly bounded 
energy E(ui,gi) < Eq < oo with respect to metrics gi = h^. ^G^ for which sup l&^l < 1. 
Then the traces Ui\gc 0 are equicontinuous. 


Proof of Corollary A.2. Pulling back the maps and metrics with the diffeomorphism 
/VV one can reduce this Corollary to the corresponding claim for the metrics Gn i and 
for maps Ui so that the functions <pf that describe the traces Ui\gc 0 ± = o ipf are such 
that there are points 9k, k = 0, 1 , 2 so that 

~ 27T ~ ~ 

(fi(9 k ) =—k and \9k+i - 9k\> C, k = 0,1,2 (A.l) 

where c > 0 depends only on 1 — sup | 6 ± | and where 6 3 := 6 q + 27 t . 

We then remark that the upper bound on i given by (3.11) implies that the metrics 
induced on the boundary of (Co,G^) are all equivalent and that the numbers Y(£) are 
bounded away from zero. Given any point p = (±1,0) we can thus apply Lemma A.l 
on a neighbourhood that is described by {(±1,0)} =F _DrO)± in collar coordinates for 
a radius r that depends only on the upper bounds on £ and l&^l, namely is chosen so 
that r < c/2, the constant of (A.l). The proof then follows by a standard argument: 
Given that the parametrisations are weakly monotone and that (A.l) does not permit 
that more than one of the three points a±(9k) is contained in the image of the small 
arc (s, 0 ) £ {± 1 } x [0 — r, 0 ± r], we then obtain the desired bound on the modulus of 
continuity from the Courant-Lebesgue Lemma. □ 

The above lemma implies in particular that any map u that is obtained as weak H 1 limit 
of a sequence of maps Ui £ Hf „(Co) is again an element of Hf ^(Cq). 

In order to pass to the limit in the differential inequality (3.3) we use at several points 
in the paper that the tangent cones Tff Hf *(Cq) depend continuously on u namely that 
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Lemma A.3. Let u € *(Cq) and let u l € *(Co) be any sequence that converges 

strongly in iP(Co) to u. Then any element v G T+ t (Co) can be approximated by 
elements in v l E T+ffp t (Co) in the sense that 

v l —> v strongly in if 1 (Co). 


Indeed, writing u\qc± = ot± o ip± respectively u % \gc± = «± ° T± and v\gc± = A± • 

a ±(‘P±) ‘ (V ; ± — < P±) we can use that strong H 1 convergence of the maps implies strong 
1 

H 2 convergence of the traces and thus also of v l \gc± '■= A± • a'±{ip±) ■ (ijj± — <p l ±) to 
v\ec±- The desired elements of *(Co) are then obtained as harmonic extensions 

of these traces similarly to the proof of Lemma 2.1 in [ 6 ]. 

As a consequence we obtain 

Corollary A.4. Let ( Ui, gi,fi ) € (Co) x M x L 2 (Cq) be such that (3.3) is satisfied 
and assume that gi —> g, Ui —» Uoo strongly in H l (Co) and fi —*■ f weakly in L 2 . Then 
(3.3) is satisfied also for the limit (u,g,f). 

We finally outline how Proposition 3.5 can be derived from the corresponding estimates 
for maps from the disc proven by Duzaar and Scheven in [ 6 , Theorem 8.3] 


Sketch of the proof of Proposition 3.5. Because of the interior estimates of Lemma 3.3 it 
is sufficient to consider points p that are contained in a neighbourhood of the boundary 
curves dC±. We then pull back the maps and metrics by a conformal diffeomorphism 
(obtained by composing with a fixed map) that maps a neighbourhood of dD\ C 
{D\ (0),g e uci) to a neighbourhood dC± C Co. As the new map u might no longer satisfy 
the three-point-condition we then modify this new triple by pulling-back with the Mobius 
transform M b ±^± to obtain a new triple (ft, g , /) = g , /) for which equation (3.3) is 

satisfied now for variations supported in a neighbourhood of the corresponding point of 
the disc. We remark that the conformal factor of g = Xg eU ci is bounded uniformly since 
we have assumed that 1 — |5 ± | is bounded away from zero and since we only consider a 
neighbourhood of the boundary where the conformal factor pi is controlled even if ^ —>■ 0 . 
Given that (3.3) holds true also for (u,g eU ci,f A 2 ) we can then apply Theorem 8.3 of [ 6 ] to 
obtain the claimed estimates on balls contained in the Euclidean disc. Since the uniform 
control on the metric Gi away from the central geodesic allows us not only to control the 
conformal factor (which appears with a different power on the left-hand side of (3.14) 
than on the right-hand side) but furthermore means that we can cover each geodesic 
ball B-j. by a fixed number of sets ip(Df(pi)) for which also the image of the 

corresponding subset of D\ with twice the radius is contained in i3f r , this implies the 
claim for the original maps. 


□ 


A.2 Properties of hyperbolic collars and the horizontal family of 
metrics Gi 

In this part of the appendix we collect some properties of hyperbolic collars, where we 
refer to the appendix of [14] and the references therein for more information, as well 
as properties of the hyperbolic cylinders (Co,Ge) that are used throughout the paper. 
We furthermore give the proof that the family of metrics described in Lemma 2.4 is 
horizontal, i.e. that fjjGe 6 Re(fH(Ge)). 
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We first recall that for S < arsinh(l) the d-thin part of the hyperbolic cylinder is described 
in collar coordinates (s,0) £ {—Y{t),Y{t)) by 


(-min(X 5 (£),y(£)),min(X 5 (£),F(£))) x S 1 


where 


X s {£) = 


2ir / n 


— arcsm 


sinh(|) 
sinh 6 


for S > t/2, respectively zero for smaller values of <5. 


(A.2) 

(A.3) 


For the metrics Gi = f^(pe(ds 2 + d,6 2 )) this means that for each S > 0 there exists a 
number co(<5) > 0 with cq(S) — > 0 for 5 —> 0 so that d-thin(Co, Ge) is contained in the 
fixed small cylinder (—co(<5), co(d)) x .S' 1 with respect to the fixed coordinates (x,6) of 
Co; or, said differently, for every c\ > 0 there exists a number S(c i) > 0 so that 

inj Q e (x,0) > <5(ci) for all \x\ > c\. 

In particular, the conformal factor p o st is bounded away from zero uniformly in i for 
\x\ > Ci. 


We also use that the norms of dz 2 on ([— Y{£),Y{£)} x S 1 , p 2 (ds 2 + d6 2 )) are given by 



II 7 2 II 87f2 

\\dz || L =o = — 

(A.4) 

and 



\\ dz IIL 2 = £3 ' 

[sin (atan(r;£)) • cos (atan(? 7 ^)) + (f — atan(r/£))]. 

(A.5) 

For £ small we thus have 

that 



\\ dz2 \\h = ^r + o(i), 

(A.6) 

while for £ large 

\\ dz2 \\h = ^ + 0(£- 5 ) 

(A.7) 


We also recall the well known fact that if a metric g evolves by dtg = Re(\ k) for a 
holomorphic quadratic differential 'k then the length of the central geodesic changes by 

dl 2n 2 

_ = -—Re(c „), (A.8) 

where codz 2 is the principal part in the Fourier expansion of or in our case simply the 
coefficient in = aodz 2 ,ao £ R. 


For large values of £, say £> Lq we can thus bound the evolution of i along a horizontal 
curve by 


M 2t r 2 \\d t g \\ L 2 

dt ~ £ ||<iz 2 ||£2 

while for small values of £ we only obtain that 


< C ■ £\\d t g \\ L 2 


(A.9) 


|^| < C ■ t^WdtgWv 

which allows for a degeneration of the metric along a curve of finite length. 


(A.10) 


Finally we explain how the formula for the horizontal families of metrics in Ai -1 claimed 
in Lemma 2.4 can be derived. 
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Proof of Lemma 2.f. Let t H > g(t ) be a curve of metrics in which moves in hor¬ 

izontal direction i.e. so that ^g(i) £ Re(H(Co, g(t))) and so that g is given as pull¬ 
back of a collar ^(— Y(£), Y(£)) x S 1 , pi(s) 2 (ds 2 + dd 2 )^ by a suitable diffeomorphism 
ft : Co —> (—Y(£),Y(£)) x S 1 where both Y(£) and fa need to be determined. 

To begin with, we derive a differential equation for Y {£) by computing the evolution of 
the width 

w{£(t)) := dist g(t) ({—1} x 5' 1 , {1} x S' 1 ) 
of the cylinder (Co,g(t)). 

Let t be any fixed time and let (s,8) £ [— Y(£(t)), Y(£{t))\ x S 1 be the corresponding 
collar coordinates. Then in these fixed coordinates, the evolution of g at time t is given 
by a o (ds 2 — dd 2 ) where ao is related to the evolution of the length £ of the central geodesic 
by (A.8). 

Thus the width of the collar, which at time t is simply given by the length of the geodesics 
m (s,Bq), evolves according to 

d d r YW , r YW 

—w(£(t)) = — / ( g ss (t)) 1/2 ds= ( g ss (t )) 1/2 • d t g S8 {t)ds 

at at J~Y(i) Jo 

r YW - i Mrf a 0 -2n f Y W , (A . U) 

= 0-0 I Pi ( S)ds = —-— / cos(— • s)ds v ’ 

Jo " Jo ^ 

= (f ) 2 «osi»(Sy(Q) = n(Ty W) |. 

For V(£) chosen so that Y(£) = — V(£)) the above formula reduces to 

d ff = -LMv W) . 


On the other hand, we can directly compute w(£(t)) by working in collar coordinates of 
g(t) as 


r Y(e) r Y(t) 

{£) = 2 / pe{s)ds = 2 / 

Jo Jo 


o 27 TCOS (A-s) 


-ds 


= 2 h(J-Y(£)) = 2/i(f - V(£)) 
where h(x) := log(tan(| + f)) is so that h'(x) = 

Thus 


dw 


= -2 


1 


di sin(V(0) dt 


(V(£)) 


meaning that V satisfies 


or equivalently 


2V' 


sin(2V) 

Thus V(£) = atan(co • £) and therefore 


= £~\ 


Z7T 7r 

Y (0 = — (2 - atan M)) 


(A.12) 
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for some constant cq > 0 . 


We can argue similarly to derive the formula for the diffeomorphism fe(x, 9) = (se(x), 9). 
Namely, we use that d t g = ff(ao(ds 2 — d9 2 )) needs to agree with 

d t9=ft' MfeiPe ■ {ds 2 + d9 2 )) = [pj{st(x)) • ((f^) 2 ' dx 2 +d9 2 )\. 


Comparing the two expressions for (dtg)ee immediately yields the condition that 

i 2n2 d i 2 \ 

and thus that for each x there exists a constant c(x) so that 


, t ... c(x) 
tan(—s e (x)) = 

Finally one can determine c(x) so that se 0 (x) = x for the number > 0 for which 
Y(to) = 1 and check that for the resulting map ft also the two expressions for ( dtg) X x 
agree. □ 


We furthermore remark that away from the central geodesic { 0 } x S 1 the metrics Gt 
converge locally smoothly as i —> 0 towards a limiting metric Go which is given by 

G 0 |c ± = r±(po(s) 2 (ds 2 + d9 2 )), (A.13) 

with f± : C± —> [0, oo) x S 1 given by 

f±(x,9) = (lim Y(£) =F s e (x),0) = • tan =F - 2?rr/,0) 

^->0 1 0 2 2 7T 

and ([0, oo), p^ds 2 + d9) 2 ) the hyperbolic cusp described in Theorem 2.7, case II. 


A.3 Properties of the diffeomorphisms h b ^ 

Here we provide (a sketch of) the proof of the properties of the diffeomorphism h b ,0 
introduced in (4.2). 


Proof of Lemma f.3. We recall that for x > 0 we can write h b ^(x,9) = (x, Xi{x)f b + + 
(1 — Ai(x))0 + A 2 {x)(j) + ) so that = A 2 (x) • is a Killing field on supp(Ai) which 
implies that 

L r lb+[ G _L L y ^G as well as Ly^G Y Ly^+G. 

To prove the other orthogonality relation we recall that since a different choice of <f> only 
results in a constant rotation on supp(YArgb) = supp(Yjb|) it is enough to consider the 
case (j> = 0 and that all vectorfields have the form Y = Y e ■ We then claim that, 
writing for short i/’ = Arg(b) and f b = f b , 0 , 

Y Arg(b)( h b(x,ip + 9)) = Yl g{b) {h b {x^-9)) while Y^(h b (x,if + 9)) = -Y^{h b {x,^-9)). 

Given that the conformal factor of the collar metric is independent of 9 this then imme¬ 
diately results in the claimed orthogonality of Ly Arg(b) G and Ly |b| G. 
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To prove the symmetry relation for Y Arg ^ or equivalently for dA ^ g ^ b ) fb we first observe 
that for b = a ■ e^, a £ M we have M b {e i ^ ,+e '>) = e^M a {e ie ) and thus 

f b ( 0 )=^ + f a ( 9-4 ,). (A.14) 


In particular 


dArg(6) 


fb(e) = l-dgf a ( 0 -ll>) 


where we can compute the derivative on the right hand side by differentiating the relation 
M a (e w ) = as 


1 _ 

1 — a 2 


d 9 f a ( 0 ) =( l M a (e lS ))- 1 • M a )(e ld ) • i ■ e l(> 


(A.15) 


(a • cosO + l) 2 + a 2 sin 2 (9 ) ’ 


Thus indeed 


fb (0 + VO = 1 - d e / a (0) = 1 - d g f a {- 9 ) = d /h(V> - 0) 


dArg(fe) ‘ 


dArg(&) 

which implies the claimed symmetry of Y Arg ( b y 
On the other hand, again for a £ R 

£fa(0) = {iM a (e ie ))~ 1 ■ (^-M a (e ie )) = - —- 2 fl ^ 2 , 2 

aa aa (1 + acos0) 2 + a 2 sin (0) 

so that for b = ae 1 ^ 

m h{ * +e)= i m= -i u - e) =-m h{ *- e) 


as claimed. 


We finally need to prove the estimates for Ly^G and Ly [fc G where we begin with the 
former for which we can use that Y rb has the simple form Y^ = A 2 (x) . 

Thus 

Ly + Ge = X' 2 (x)p 2 (sg(x)) ■ (dx <g> d6 + d9 <g) dx) 

and as pg o sg is bounded uniformly both from above and below on any fixed cylinder 
[d, 1] x S 1 C Go for l £ (0, Lo] we easily obtain the claimed estimate for Ly^G. 

To analyse Ly |b| G we first remark that 

K,<p G e = Pe(s e (x)) ■ [(^) 2 + (% rY] dx 2 + (%*) ■ {~cw -) • ( dx ® dd + d9 ® dx) 

+ (^) 2 ^ 2 ] 

(A.16) 

and that in view of (A. 14) we only need to consider the case that b = a £ R. 

As we only wish to prove a lower bound on \\Ly w G||x ,2 it is enough to consider the 
subcylinder [|, 1] x S 1 on which the above expression reduces to 

KGg = p 2 (sg(x)) ■ [(f|f fdx 2 + (^) 2 d0 2 ] 
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so that 


±(h* a G) = 2p 2 (se(x))-^-£^de 2 . 

On this part of the cylinder we furthermore have that g ee = p ~ 2 o sg • f )~ 2 where we 

recall that p is again bounded uniformly from below so that 


\Ly w G \\Uco,g) = \\£(KG)\\l*(c 0 ,h*g) >4 

>c • 


'?fi\- 2 \jL0U\ 2 dVo 


[|'l] xS 1 


1£4V 
v de ) 


' da dO 


f (£ik) 1 \ 2 d0 

/ V de ) \da de \ aa 

Js 1 

1 — 2 

for some fixed constant c > 0. Now as dgf a = (i +acos gy° +Q a sin i^ g we can compute 
dadefa = ~ \ae lff + 1| 4 - 2a(l + acosf?) 2 + 2a 3 sin 2 9 

+ (1 — a 2 ) [2 cos(0)(l + a cos 9) + 2a sin 2 0] 

= — | ae lB + 1| 4 • [2a sin 2 6 + 2(1 + a cos 6) ■ (a + cos#)] 


(A.17) 


(A.18) 


We set £ = 1 — a and remark that for e small and for 9 given by 9 = n + A ■ e 
[2a sin 2 9 + 2(1 + a cos 9) ■ (a + cos 0)] > 2 • [A 2 £ 2 — £ 2 + 0(e 3 )] 
so that it is in particular bounded away from 0 by 2s 2 for angles 2e < \9 + 7r| < 3 s. 
Combined with (A.17) we thus find that 

r n+3e 1 _ a 2 


\\ L Y m G\\ L 2( CojG) >C- £ 


& 


ae l9 + 11 “ d9 > ce 2 = 


> 7r+2e |ae* e + 1| 

for 1 — a sufficiently small. This implies the claim of Lemma 4.3. 


(1 — a) 2 

(A.19) 


□ 


Sketch of Proof of Lemma J^.2. The property asked for in Lemma 4.2 is esstentially a 
consequence of us choosing the diffeomorphisms as restrictions of Mobius transforms 
onto S 1 and the fact that given any two triples (wi,W 2 ,W 3 ) and (z\,Z 2 ,Z 3 ) of points 
on S 1 there is a unique Mobiustransform mapping Zi to Wi. To be more precise, using 
the group property of the Mobius transforms one can reduce the claim of Lemma 4.2 to 
proving that for any distinct ^ 1 , 2,3 £ [ 0 , 27r) and any a 0 £ [ 0 , 1 ) the derivative of the map 
(h,%f) (/b^('di), fb^i^), fb^i^)) has full rank in the point {b,ip) = (a o ,0) £ C x R. 

A short calculation then verifies this claim. □ 
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